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Abstract 

An n-body system is a labelled collection of n point masses in a Eu- 
clidean space, and their congruence and internal symmetry properties in- 
volve a rich mathematical structure which is investigated in the frame- 
work of equivariant Riemannian geometry. Some basic concepts are n- 
configuration, configuration space, internal space, shape space, Jacobi 
transformation and weighted root system. The latter is a generalization 
of the root system of SU(n), which provides a bookkeeping for expressing 
the mutual distances of the point masses in terms of the Jacobi vectors. 
Moreover, its application to the study of collinear central n-configurations 
yields a simple proof of Moulton's enumeration formula. A major topic is 
the study of matrix spaces representing the shape space of n-body configu- 
rations in Euclidean k-space, the structure of the m-universal shape space 
and its 0(m)-equivariant linear model. This also leads to those "orbital 
fibrations" where SO(m) or 0(m) act on a sphere with a sphere as orbit 
space. A few of these examples are encountered in the literature, e.g. the 
special case /0{2) k, was analyzed independently by Arnold, Kuiper 
and Massey in the 1970's. 
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1 Introduction 

An n-body system is defined to be a labelled collection of n point masses (or 
particles) Pi of mass rrii in Euclidean 3-space M^, and it is of general importance 
to find appropriate mathematical models to describe and analyze such a sys- 
tem. We have in mind few-body systems as well as many-body systems, ranging 
from different areas such as celestial and quantum mechanics or quantum chem- 
istry (n-body problem, n-atomic molecules etc.). Despite the diversity of the 
applications they all share a fundamental underlying mathematical structure, 
in terms of kinematic concepts and internal space geometry, and the role of 
the mass distribution. Here we shall focus attention on these basic structures, 
in a modern geometric and topological setting with orthogonal transformation 
groups in the forefront. This approach also establishes similar results for point 
masses in any Euclidean spaces M"', d > 2. 

In this introductory section we first give an overview of the paper, which 
take up several different topics. Then, in the following subsections we introduce 
some basic concepts and constructions which we shall return to later. Here 
the presentation is rather informal or expository, with some comments on the 
history. 
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1.1 A brief overview 

Our approach is to combine symmetry and kinematic geometric principles in the 
framework of equivariant geometry, exhibiting the importance of the classical 
orthogonal transformation groups and the associated orbit space constructions. 
This enables us to investigate in a unifying way the notions of congruence, inter- 
nal configuration space and internal symmetry group. Internal symmetries of 
n-body systems are investigated in Section 2 from a categorical viewpoint based 
on pure geometrical principles. As we shall explain, for various reasons they 
should be referred to as Jacobi transformations, but in the physics literature 
they are also encountered as kinematic rotations or democracy transformations, 

cf. nzi. 

Section 3 is largely devoted to the study of the topology of shape spaces in 
general, and here it is natural to consider n-body systems in higher dimensional 
Euclidean spaces as well. We shall exploit the fact that congruence and 
internal symmetry for point masses in M.'^ combine together to a tensor product 
representation of some 0{d) x 0{m), acting on the matrix space M{d, m) ~ W^"^ 
by matrix multiplication. In short, a typical shape space is the orbit space 
M* = SP/0{d), where is the unit sphere of M{d,m). Now M* inherits the 
action of 0{m) as its symmetry group, but on the other hand, M* also naturally 
embeds as an 0(m)-invariant subset of the linear space Sym?{m) of symmetric 
matrices with zero trace. Of particular interest are the cases where S^/H ^ S*^ 
is also a sphere, H = SO{m) or 0{m). They yield an infinite family of "orbital 
fibrations" S'^, which for q = 1,2,4 are the Hopf fibrations, cf. Section 

3.4.2. 

In the study of many body systems, we apply the orbit space reduction to 
the configuration manifold rather than to its cotangent bundle, as in the re- 
duction method of Marsden-Weinstein which, for example, Iwai ^H] applies to 
the Hamiltonian system describing classical molecular dynamics. On the other 
hand, R. Littlejohn and his collaborators (cf. e.g. jl6j. jl7j . jl8j ) have investi- 
gated the gauge fields that arise on the reduced (i.e. internal) configuration 
space, and our work is directly related to the geometric framework of their 
investigations. 

In Section 4 we determine the geometric invariants of n-body systems, 
namely polynomial functions on the configuration space which are both con- 
gruence invariants and internal symmetry (or democracy) invariants. The topic 
is certainly well understood in classical invariant theory, but we are also seek- 
ing symmetrized expressions for centered n-body configurations, that is, their 
center of mass is fixed at the origin. 

In Section 5 we introduce the notion of the weighted root system of an 
n-body system with a given mass distribution. This is a geometric invariant 
which generalizes the notion of a root system of Cartan type An-i in classical 
Lie theory, and its underlying structure is actually inherent in various contexts. 
For example, it encodes the data of the relative positions of the binary collision 
varieties in the configuration space or shape space. We shall introduce it as a 
bookkeeping device for expressing the mutual distances between the n bodies 
in terms of the Jacobi vectors. 
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In Section 6 we present, as a nice application of the weighted root system, 
a simple proof of Moulton's classification of collinear central n- configurations 
(cf. [23]) which is also well adapted for numerical computations. Recall that 
these are the configurations characterizing n collinear masses capable of a rigid 
uniform rotation under the mutual gravitational forces. 

The present paper is essentially a preprint (with the same title) of the au- 
thor from April 2002, prompted by the paper Hsiang and the succeeding 
joint work Our Sections 5-6 recall the contents of but in the present 
paper the topology of general shape spaces, the universal shape space and its 
linear model, is a major topic. It turns out that some few-body shape spaces 
have an interesting history in the literature. For example, it was surprising to 
find that the quotient space of complex projective plane CP^ modulo complex 
conjugation is topologically a 4-sphere, namely 5^/0(2) S^. Several indepen- 
dent and different proofs of this fact had already been published (cf. Arnold 
[2], [SI, Kuiper[T^, Masseyj^^). Now we also find that our Section 3 has some 
overlapping with the more recent paper Atiyah-Berndt [T|. 

1.2 n-configuration space, internal space and shape space 

The location of an n-body system is conveniently represented by its n-configuration 
X = (a^,a2, ..,a„) where is the position vector of Pi. The n-tuple X is re- 
garded as a vector in the free n-configuration space 

Mn = X X .... X ~ IR3n^ ^^-^ 

namely a Euclidean 3n-space with an orthogonal splitting reflecting the in- 
dividual positions of the n point masses, which may possibly coincide. The 
mass distribution (mi, m2, •., rra„) is tacitly assumed to be fixed unless other- 
wise stated, and in general we assume mj > 0. The centered n-configuration 
space is the subspace of dimension 3n — 3 

consisting of those n-configurations with its center of mass at the origin. Mn 
has the following mass dependent (Jacobi) kinematic metric, namely the inner 
product of X and Y = (b-^, b2, .., b„) is 

(X,Y) = ^m,a,-b, (3) 

The isometry group of Mn is the associated Euclidean group 

E{Mn) = M„,xO(M„) ~ R^"xO(3n) = E{3n) (4) 

which we have expressed in the usual way as a semidirect product of the sub- 
groups of translations and orthogonal transformations respectively. 

Two n-configurations X and Y are regarded as congruent if their i-th com- 
ponents a, and bj differ by the same Euclidean motion in for each i. Thus 
the congruence relation is defined by the natural (diagonal) i?(3)-action on M„, 
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by which E{3) (resp. 0(3)) is embedded as a subgroup of E{3n) (resp. 0(3n)). 
The space Mn of congruence classes is the (congruence) moduli space, also re- 
ferred to as the internal space. It may be viewed as the result of the two-step 
orbit space constructionb 

where in the first step the translation-reduced space is identified with M„. 
The second step is the orbit space construction of the transformation group 
(0(3), M„), from which it follows that M„ is a stratified manifold of dimension 
3n — 6 (cf. Section 3). 

Next, let us divide "congruence" into "shape" and "size" and use the squared 
norm function / (cf. ©), namely the polar moment of inertia, as a measure 
of the size of an n-configuration (or of its congruence class). Then M„ has the 
structure of a cone, where each ray emanates from the cone vertex (or base 
point) O = (/ = 0) and the ray represents a fixed shape. The shape space is 
the space of rays, which we may regard as the orbit space 

Mn-{0} M„-{0} ^ 

M+ 0(3) X M+ " " ^ ' 

where 0(3) x M"*" is the group of similarity transformations of Euclidean 3-space 
with the induced (diagonal) action on M„. However, each ray has a unique point 
where 1 = 1, and therefore it is more convenient to identify M* with the subset 
of Mn consisting of classes of unit size 7 = 1. Thus, the internal space M„ is 
naturally a cone over the shape space M*. 

1.3 Geometrization, symmetry and reduction 

The viewpoint that kinematics is a geometric discipline has a long history. For 
an n-body motion t X(t), the fundamental kinematic quantities are 



|2 I |2 

mi a,; 



2 



X| 

n = X X X 



m,a. 



These are the moment of inertia, kinetic energy, angular and linear momentum, 
respectively. In particular, the hyperradius p = \fl is the norm of X in the 
n-configuration space (Q) with respect to the inner product Q, and using T 
the same metric may be presented as the kinematic Riemannian metric 

ds^ = 2Tdt^ = ^ nii{dxf + dyf + dzf) (8) 

Furthermore, dynamics was incorporated in this geometric setting by the 
classical geometrization procedure which dates at least back to the early 19th 
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century. Let us briefly recall the basic idea behind this, namely Jacobi's re- 
formulation of Lagrange's least action principle, which goes as follows. When 
dynamics is taken into account and the above n-body motion is due to a force 
field VC/ derived from the potential energy —U, the total energy h = T — U is 
conserved and the metric © should be conformally modified to the following 
dynamical metric 

dsl = {U + h)ds^ (9) 

depending on a fixed energy level h. Then the trajectories of total energy h 
can be recovered as the geodesies of this metric. This applies, for example, 
to physically important n-body systems with potential functions (such as the 
Newtonian or Coulomb potential) depending on the pairwise distances rjj = 
|aj — SLj\. Therefore, they are invariant under the Euclidean group E(3), and 
consequently both linear and angular momentum of the motion are conserved. 

Since the early days it has been an important issue how to fully utilize 
conservation laws. Sometimes these are first integrals associated with sym- 
metry groups, and the reduction of integration problems using continuous (or 
infinitesimal) symmetries dates back to Sophus Lie's work in the 1870's. 

Clearly, the most effective usage of the invariance of linear momentum p 
is to choose an inertial frame at the center of mass, thereby reducing the con- 
figuration space M„ to its subspace M„ ((2} and hence the associated linear 
momentum vanishes. A much harder problem is to further use the invariance 
of ft, which is related to the congruence action of 0(3) (or 50(3)) on M„ and 
its cotangent bundle. Namely, if T* denotes the cotangent bundle construction, 
consider the two ways of reducing the phase space 

^ ^0{3)'~ " ' 0(3) 

to dimension 6n — 12 and 6n — 9, respectively. This suggests that the orbit 
space reduction is most effective at the level of the configuration space rather 
than the phase space. 

In this paper we shall focus attention on two main issues. Firstly, it is im- 
portant to construct appropriate coodinates for M„, which is actually a problem 
with no canonical or generally "best" solution. However, from the representa- 
tion theory of 0(3) and natural guidelines suggested by the splitting or invari- 
ance of the kinematic quantities, we are led to a natural approach whose origin 
may, in fact, be ascribed to Jacobi. Secondly, we inquire about the structure of 
the internal space Mn, as an orbit space of M„ modulo 0(3). It has the natural 
kinematic metric 

df = 2f dt^ (10) 

where T = T — T"^ and T"^ is the kinetic energy due to purely rotational 
motion. Indeed, this metric is Riemannian and coincides with the induced 
orbital distance metric (see e.g. jSO])) and the reduction map M„ Mn is a 
(stratified) Riemannian submersion. As a consequence, the geodesies in M„ are 
the image curves of those geodesies (i.e., linear motions of constant speed) in 
{Mn^ds^) with vanishing angular momentum. 
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On the other hand, the appropriate reduction of the dynamical equations 
in Mn to the level of M„ is generally a hard problem even today. In M„ the 
geometrization procedure yields the conformal modification of the kinematic 
metric © leading to the dynamical metric ©. Similarly, one may search for a 
similar geometrization procedure at the internal space level, which should yield 
a dynamical Riemannian metric on M„ of type 

dsln = F{U,h,Q)ds^ (11) 

in analogy with Q. However, this only works in special cases, for example, for 
the (classical) planar n-body problem, but we leave this topic here. 

Actually, in his study of celestial mechanics Jacobi himself abandoned the 
above Riemannian geometric approach in favor of the increasingly successful 
Hamiltonian formalism where, at the phase space level (symplectic geometry), 
ideas involving symmetry, conservation laws and integrability questions have 
been continuously developed up to present time. On the other hand, with 
the modern techniques of equivariant geometry, Lie transformation groups and 
related reduction theory and quotient constructions, the framework of Rie- 
mannian geometry and modern differential geometric techniques are now more 
readily applicable for the study of n-body problems. We also refer to Littlejohn- 
Reinsch[T7) for a general discussion with many references to the gauge kinematic 
and dynamics of many particle systems. For real historical background infor- 
mation on the geometrization of physics we also propose the survey article 
Liitzen^^ on the 19th century interactions between mechanics and differential 
geometry. 

1.4 Congruence and internal symmetry in n-body spaces 

Let us explain the interaction of the notions of congruence and internal sym- 
metry for n-body systems, where by congruence we mean 0(3)-congruence and 
assume the translational degrees of freedom have been eliminated. For this pur- 
pose we introduce a notion somewhat more general than a free m-configuration 
space (P), namely an m-body space 

M-M^m . (o(3),mp3) 

is a Euclidean 3m-dimensional space with an orthogonal representation of type 
(0(3), mp^), that is, m copies of the standard representation P3 of 0(3). How- 
ever, particles and mass distributions are not mentioned, and there is no spe- 
cific decomposition into 3-dimensional invariant subspaces M.^ as in As an 
example, we have in mind the centered n-configuration space M„ 0, where 
m = n — 1 and the original mass distribution is disguised in the Euclidean 
metric ©. 

On the other hand, associated with M^™ is the totality P(M^™) of all or- 
thogonal, 0(3)-invariant and ordered decompositions 

m3™ = e ^2 e . . .®Vm , Vic^ R^. (12) 

Then, for a specific choice of decomposition we can identify R^™ with the free 
m-configuration space Mm as in and associate a particle Pi with position 
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vector aj to each summand Vi. Moreover, the metric on M^"* viewed as a metric 
on Mm should have the form so we also need a mass distribution. However, 
by suitably scaling of the vectors aj we are actually free to choose any mass 
distribution since a vector in M'^™ written as an m-configuration (bi, .., b^) 
determines position vector aj and mass mj modulo the constraint bj = y^mjaj 
for each i (see Section 2.1). 

In any case, R^™' has the given conguence group 0(3) C 0(3m), and we de- 
fine the internal (or inner) symmetries to be those transformations ip £ 0{3m) 
which commute with 0(3), that is, ipip = iptp for each £ 0(3). By standard 
representation theory (cf. Schur's lemma) these ip constitute a subgroup iso- 
morphic to 0{m), which we shall refer to as the (internal) symmetry group. 
The two subgroups of 0{3m) intersect at Z2 = {it/d} and hence combine to a 
subgroup 

^ 0(3m) 

£12 

which is also described mathematically by the tensor product p^ of their 
standard representations. 

We also point out that 0{m) acts naturally on the set of decompositions 
(|12|) . This action is, in fact, transitive and has isotropy group of type O(l)™', 
and hence establishes a 1-1-correspondence 



I?(M-^'") 



s ^, 0{m) 



0(1) 



To be more explicit, let us fix some decomposition (|12|) and choose orthonor- 
mal bases to identify each Vi with M^, and hence identify R^'" in H12|l with the 
space M{3, m) of real 3 x m- matrices 

X = (xi,X2,..,Xm) (13) 
with the standard Euclidean norm square 

m 

\X\^ = trace{X^X) = ^\^i\^ , (14) 

i=l 

where Xj G Vj is the i-th column vector of X. Then the action oi tp £ 0(3) and 
if £ 0{m) is just matrix multiplication on the left and right side respectively, 
inducing a joint left action on matrices by 

{'il).,Lp)X = iI}Xlp^'^ (matrix multiplication) (15) 

Thus congruence and symmetry combined together is the following tensor prod- 
uct representation 

(0(3) X 0(m),P3®p„,M(3,m)), (16) 

whose orbit structure will be analyzed by combining two consecutive orbit space 
constructions 

M(3,m) , M(3,m) M(3,m) , 
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Of primary interest is the first orbit space, namely the internal space M = 
M(3,m) with its induced metric structure. The symmetry group 0{m) de- 
scends faithfuhy to an induced transformation group on M which is, in fact, its 
isometry group. Thus, the last step in H17() yields the final orbit space which is 
geometrically the Riemannian cone over a spherical triangle A (resp. a circular 
arc if m = 2) : 



Of course, this is a Euclidean cone. However, as indicated it is also the funda- 
mental domain of the Weyl group B3 (resp. B2 ii m = 2) in classical Lie theory, 
and this tells us that C(A) embeds isometrically into M(3, m) as a cross section 
of the transformation group (|16|) . Namely, C(A) hits every orbit of that group 
at a unique point and, moreover, it is perpendicular to every orbit. 

We refer to Section 3 for further analysis of the above orbit spaces, in a more 
general setting involving all matrix spaces M{d, m) with the transformation 
groups 0{d) X 0{m), for any d> 1. 

1.5 Jacobi vectors and the centered configuration space 

In dynamics it is the conservation of linear momentum p that enables one to 
reduce the n-body problem to an (n-l)-body problem plus a trivial 1-body prob- 
lem for the motion of the center of mass. To explain this, consider the canonical 
orthogonal and 0(3)-invariant decomposition of the free n-configuration space 



where the subspace M„ ~ IR'^"^^ is the centered n-configuration space ((SJ and 
the "diagonal" AM'^ = {(a, a, .., a)} is its orthogonal complement. For a motion 
X(t) in Mfi the vector a = a(t) is, indeed, the center of mass of X(t). Now, 
conservation of p means our inertial frame of reference will remain inertial if we 
translate the frame so that a becomes the new origin. Therefore, with respect 
to the new frame, the motion will take place in the summand M„ in 
This simple reduction is a key step in the integration of the classical Kepler 
problem, where n = 2 and M2 ~ ~ Mi is the configuration space of a 
ficticious 1-body system. 

However, for n > 2 M„ is just an (n-l)-body space and there is no canonical 
way of further decomposing M„ to become the configuration space of n — 1 
ficticious particles with appropriate position vectors and mass distribution, cf. 
1)12(1 . This fact is clearly reflected by the variety of types of coordinates for 
M„ which can be found in the literature. Recall, for example, the efforts of 
Lagrange, Jacobi and Delaunay who constructed their own "good" coordinates 
to study the classical 3-body problem (cf. e.g. Marchal j20l|). 

Of particular interest to us is Jacobi's approach, which by repeated appli- 
cations generalizes to n bodies, but again there is no canonical way of doing 
so. Anyhow, our interpretation of his basic idea is that a solution of the above 




(18) 



M„ = M„ e AM^ , cf. JH) 



(19) 
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splitting problem amounts to the construction of a transformation 

^: Mn^Mn-l (20) 
: (a-|^, a2, .., a„) — > (xi, X2, .., x„_]^) 

with the "appropriate" properties (see below), connecting the (n-l)-body space 
M„ (with the Jacobi metric ©) to a standard model, namely the free configu- 
ration space Mn-i with all masses equal to 1. 

By viewing the vectors Xj G as the columns of a matrix we may 
identify Mn-i with the matrix space 

n-l 3 

M(3, n - 1) = M"-^ = ^ eMf = ^ eM^-^ (21) 

i=i 3=1 

with the norm as in (|14() . congruence group 0(3), symmetry group 0{n — 1) 
and their joint tensor product representation of 0(3) x 0(n — 1), see (|15p. (jl6p . 
In H21|) we have also indicated that the column and row vectors of a matrix 
belong to two different Euclidean spaces, having the standard action of 0(3) 
and 0{n — 1) respectively. 

Now, what are those "appropriate" properties Jacobi transformations such 
as ^ should satisfy ? In short, the answer is that ^ is just an 0(3)-equivariant 
isometry, that is, 

i) ^{gX) = 5^(X), g G 0(3), ii) |X| = |vI/(X)| (22) 

Let ^ in ()2U() be a given transformation of this type. It associates to a centered 
n-configuration X = (a]^,a2, ..,a„) its Jacobi vector matrix 

X = ^(X) = (xi,X2,..,x„_i) G M(3,n - 1) 

whose column vectors Xj will be referred to as the corresponding Jacobi vectors. 
On the other hand, if ^' is another Jacobi transformation, then the composition 
o ^I*"^ = ip is still an 0(3)-equivariant isometry. Consequently, we can write 
vj/' = o vj/ where 

V9 : M(3,n- 1) ^ M(3,n- 1) (23) 

is a Jacobi transformation of the standard model (|21|). namely an orthogonal 
transformation commuting with 0(3). In other words, if belongs to the sym- 
metry group 0{n— 1) (acting on the matrices X by right multiplication). This 
explains the non- uniqueness of Briefly, by knowing only one of the ^'s we 
obtain all of them by composing ^ with any ip G 0{n — 1). 

So far, however, we have not constructed a single ^ in (|2()|) . but we have 
just seen that this is all we need to do. We refer to Section 2.2 for the ex- 
plicit construction of our standard choice ^' = ^'o of Jacobi transformation 
and hence also our preferred choice of Jacobi vectors Xj as distinguished linear 
combinations of the vectors aj. 

Remark 1.1 Quantities which are invariant under the transformation group 
{0{n — 1),M„) are also insensitive to the different orderings of the n bodies. 
Therefore, this group is sometimes referred to as the democracy group. Ele- 
ments (f G SO{n — 1) are also named kinematic rotations. See e.g. \n\ . 
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2 Jacobi transformations from a categorical view- 
point 

Loosely speaking, there is the category of m-body spaces whose objects are 
the Euchdean spaces M ~ M^*" with an orthogonal transformation group of 
type (0(3),m/33), and the morphisms will be called Jacobi transformations, 
namely they are the 0(3)-equivariant isometrics M — > M' between the spaces. 
Here the matrix space M(3, m) = Mm is a distinguished object which also 
serves as our standard model (|21l) . Therefore, for a given m-body space M a 
"good" coordinate system for M amounts to the choice of an appropriate Jacobi 
transformation 

^ -.M ^ M{3,m) 

and the columns of the matrix X G M(3, m) will be the Jacobi vectors with 
respect to ^. We are primarily interested in the centered n-configuration space 
M = Mn (where m = n — 1). 

In the following two subsections we shall characterize Jacobi transformations 
in several equivalent ways, and we start with the classical Jacobi construction 
which also justifies our usage of terminology. 



2.1 Transformations of n-body systems and Jacobi's approach 

In the Introduction we actually started (less formally) with the category of n- 
body systems, whose objects are n-tuples (Pi,P2, ■■■,Pn) of point masses Pi = 
(aj, nii) with a, and rui as position vector and mass respectively. In this setting, 
what should be the appropriate morphisms 

<I>:(Pi,...,P„)^(P{,...,P^) ? (24) 

After all, the usefulness of such transformations with the desired properties, to 
simplify the further analysis, is well documented in both classical and quantum 
mechanics. First of all, we propose to consider transformations 

$: (ai,--,a„)^(a;....,a') 
(mi, ...,m„) ^ (mi, ...,m„) 



where the first map in (|25j) is an invertible linear transformation on M„ = M^"' 
which may possibly depend on the masses m^, with the choice of "new" masses 
constrained in some way. But $ also transforms the motions of n-body 
system and hence also the first three basic kinematic quantities in ((Jj), so let us 
demand that $ preserves them, namely 



I = ^^'TT'j |ajp = ^^w-'o la' 



2 / I / |2 

i Pi I 



2r = ^mi|ai|^ = ^m^- |a^|^ (26) 

ft = rriiai x a.j = ^'j^'j ^ ^'j 

Motivated by an idea of Jacobi, we define the (i,j)-basic Jacobi transforma- 
tion for i ^ j 

:(Pi,...,P„)^(Pi',...,P;) (27) 



11 



by demanding = Pk,m'i^ = ruk for k ^ and 



Pi = (a'„ m'J = (a, - a„ ' ' ) (28) 

rrii + rrij 

p. = (a.,m- = ( ■ '-^,mi+mj) 

■' J J rrii + nij 

In particular, Pj is the center of mass of Pj and Pj, and m'- is their reduced 
mass 

P^ij = \_ ' 29 
■' rrii + rrij 

Among the basic Jacobi transformations let us also include the following mass 
normalizing transformation 

Pq : (aj,mj) {^/m'i^^, 1), Vz (30) 

and the permutation transformations 

a: (Pi,...,P„) ^ (P<,(i),...,P,(„)) (31) 

which permute the vectors aj and masses mj covariantly, namely there is one 
for each permutation a of {1, 2, .., n}. It is straightforward to check that the 
above transformations, and hence all their compositions, have the invariance 
property ((^ . 

Next, we inquire whether the above special transformations already generate 
all those $ with the property H26() . This is affirmatively settled in the following 
subsection, but first we shall make some clarifying observations. Note that 
velocity vectors aj of a motion are transformed in the same way as the position 
vectors a^ and may well be regarded as vectors bj independent of the a^'s. In 
particular, invariance of / means the same as invariance of T, namely that $ is 
an isometry. Therefore, the remaining issue is the nature of the invariance of 
ft. 

Consider the 3n x 3n-matrix of a linear transformation <I>, regarded as an 
n X n-matrix 

(Bu Bi2 \ 
B21 B22 (32) 

where each block Bij is a 3 x 3-matrix. Such an $ expresses each a'j as a linear 
combination of the old ones, that is, 

n 

= ^ f^ijAj, i = 1, .., n, (33) 
i=i 

if and only if each 3-block is a scaling matrix Bij = f3ijld, where Id denotes 
the identity matrix. Equivalently, $ commutes with the (diagonal) action of 
GL{3) on M^" and hence (by representation theory) belongs to the subgroup 
GL[n) embedded into GL(3n) by the tensor product action of GL{3) x GL{n). 
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On the other hand, assuming ((321) holds $ will be an isometry (i.e. / is 
preserved) if and only if the n x n-matrix {Pij) is "orthogonal" in the sense that 

iPij) G D'-^0{n)D C GL{n) (34) 

where D = diag{^/rni, ^/rn^), D' = diag{y^ m'l, .., ^Jra\^. Finally, assuming 
and it is easy to verify that is also preserved. 

2.2 Equivalent characterizations of Jacobi transformations 

It turns out that the invariance properties H26|l are satisfied solely by demand- 
ing invariance of the vector f2. This is elucidated by the following proposition 
where the angular momentum construction is analyzed from a purely algebraic 
viewpoint. 

Proposition 2.1 Consider two n-tuples of vectors in 3-space 

X = (xi,x2, ..,x„),y = (yi,y2, ..,y„) 
and define their angular momentum vector by 

n 

n = n(x,y) = ^x,xy, 

i=l 

Regard X, Y as vectors in M^" = X .... X and consider general linear 
transformations <I> : M^" — > M'^", 

X ^ X' = Y ^Y' = (y'i,y2,..,yn) 

Then <I> preserves angular momentum, that is, for all X, Y 

fl{X,Y) = fl{X',Y') 

if and only if ^ is an 0{3)-equivariant isometry. 

Proof. The 3 x 3-matrices 

Sk = {sif), k= 1,2,3, 

whose only non-zero entries are s[j^ = —s^^^ = 1 whenever {i,j,k) is an even 
permutation of (1, 2, 3), constitute a basis for the Lie algebra so{3) of all skew- 
symmetric matrices 5. The action of 50(3) on M^" embeds each 5 as a "diag- 
onal block" matrix 

A(5) G so{3n) 

with n copies of S along the diagonal. Observe that the matrices A{Sk) actually 
represent the three components of the vector ft , as skew-symmetric bilinear 
forms on R^" expressed as matrix products, namely 

nk{X,Y) = X'A{Sk)Y 
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where X, Y are regarded as column vectors with 3n components. 

Let us identify a transformation $ with its matrix in GL{3n). Those trans- 
formations leaving invariant constitute a matrix group Gk defined by the 
constraint 

^'A{Sk)<^> = A{Sk) 

and hence the group leaving ft invariant is the intersection G = Gi f] G2 H G3 
defined by the constraint 

$*A(S)$ = A{S), for all S G so(3) (35) 

The proof amounts to show G = 0{n) C 0{3n), namely that G is the group of 
isometries commuting with 0(3) (or equivalently with SO (3)). 

If we knew $ is orthogonal, then we could have applied the exponential 
function to (|35() to conclude that $ commutes with 50(3) and hence belongs 
to 0(n). We claim, however, the isometry condition is itself a consequence of 
(|35)) . To see this, write a typical matrix of the Lie algebra L{G) as an n x n- 
matrix (Aij) with blocks Aij, see H32() . Then the infinitesimal version of H35() 
simply reads 

AljS + SAij = 0, for all i,j and all 5 G so{3), 

or equivalently 

Aij = —Aji, Aij = aij Id (scaling matrix) 

However, these are precisely the conditions defining the Lie subalgebra so{n) 
and consequently G and 0(n) have the same connected component SO{n). 

Finally, to show G = 0{n) we use the fact that G lies in the normalizer of 
SO{n) in GL{3n), namely in the image of GL{3) x 0{n) by the tensor product 
representation. Therefore, we may assume <I> = A(i?) where B G GL{3) and 
hence by 

B^SB = S for all S G so(3) 

This condition says B leaves invariant all skew-symmetric bilinear forms in M^, 
and consequently B = ztld. This completes the proof. ■ 

The basic Jacobi transformations (|7f)) . (PU)) . (|ST|) are maps between n- 
configuration spaces with possibly different mass distribution. However, using 
normalizing transformations Pq, as in (|3()j) . any composition may be "pulled 
back" to a transformation M(3,n) M(3,n) and with all masses rrii = 1. 
In this interpretation, the (i,j)-basic Jacobi transformation pij will be a (2- 
dimensional) rotation Pij G 0(n) and it is easy to check that Pu, P23t P13 S^^" 
erate a dense subgroup of SO (3) C 50(3) x SO{n — 3) C 0(n). Therefore, all 
rotations Pij together with the permutations (|^T|) generate a dense subgroup of 
0(n). 

As a consequence of the above observations we now state the following 
result on the characterization of Jacobi transformations, a terminology justified 
by (iii) below. In view of Section 1.3 we shall not attempt to define the angular 
momentum vector ft in the broad category of n-body spaces, but in the following 
theorem we refer to n-configuration spaces as defined in Section 1.1. 
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Theorem 2.2 The following four classes of linear transformations $ between 
n- configuration spaces are identical : 

(i) ^ is an isometry which preserves the angular momentum vector Q,. 

(a) ^ preserves fl . 

(Hi) ^ is in the closure of the set of transformations generated by the basic 
Jacobi transformations. 

(iv) ^ is an 0{3)-equivariant isometry. 

The Jacobi transformations on a fixed (n-body or n-configuration) space are 
also its symmetries (cf. Section 1.3), and they constitute a group isomorphic to 
0(n). In particular, the symmetry group for our standard model, namely the 
matrix space M(3, n), is the group 0{n) acting by matrix multiplication on the 
right side. 

Transformations which preserve angular momentum and kinetic energy are 
certainly useful in quantum mechanics. Here one tries to keep operators "sepa- 
rable", that is, with no cross terms (cf. e.g. §10.1). As a consequence of the 
above results Jacobi transformations provide all possible linear combinations a^- 
of the vectors aj which preserve the kinetic energy operator, in the sense that 



2 



±vl + ... + —vl 



2 



m 



V^a' + - + 



(36) 



A reduction of variables is achieved by choosing the transformation <I> so that 



4 



mass coordinate system. Then, for example, the (time dependent) Schrodinger 
equation with potential energy V reduces to 



5 T 



2 



n— 1 ^ 

(E— ^a') + na'i,...,a; 



n—li 



where ^ = ^'(a'^, a'^_-|^, t). In the following subsection an explicit construc- 
tion of such coordinates is given, and by usage of the symmetry group 0{n) 
they can be modified to satisfy additional properties for a specific purpose. 

On the other hand, in the study of atomic structures one of the masses, say 
m„ (the nucleus), may be relatively large and then non-Jacobian coordinates 
may also turn out to be useful. However, in that case mass polarization terms 
cannot be avoided, namely mixed terms V^' • V^' will appear in the kinetic 
energy operator (see e.g. Appendix 8 in ^ 

2.3 Construction of Jacobi vectors 



Here we shall construct an explicit Jacobi transformation (|2n|) . or rather its 
inverse 



^ ri—l 

$ = -Jf"^ : X = (xi,X2,...,x„„i) (a^,a2, ..,a„„i, V'l 



rriiai 



(37) 



i=l 
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from M(3,n — 1) to the centered configurations space M„. By definition, $ 
must be an 0(3)-equivariant isometry, and hence in analogy with H22|) 

i) ^gX)=g^X), geO{3) (38) 
ii) \^{X)\ = \X\, X £ M{3,n-l) 

and let us also add the following "normalizing" condition : 

(xi,X2,..,Xfc,0,...,0) ^ (ai,a2,..,afc,b(^'),..,b('^)), for A; < n - 1 (39) 

Namely, the Jacobi vectors Xj should vanish for i > k ii and only if the particles 
Pi "collide" for i > k, that is, they occupy the same position b'-'^-*. Equivalently, 
the Jacobi vectors define an orthogonal decomposition 

Mn = Rl eMl e ... 

where ~ is the image of the imbedding 

(0, .., 0, Xfc, 0, .., 0) ^ (0, .., 0, - b('=-^), b(fc), .., b(,)) (40) 

which actually represents a collinear 3-body system in the sense that Pi,P2, Pk-i 
are located at the origin and Pk+i^...., Pn at the common position vector 

b(,) = bW-b('=-i). 

The condition that X = (ai, .., a^, b^''), b^'')) in ((HHI) belongs to M„ yields 
the formula 

k 

b^'^ = ;^E"^*^- bW=0, (41) 

i=l 

where by definition 

m^^^ = rrik+i + ... + rUn, fn^^^ = mi = fh 

As a consequence of H39() we also note that the matrix [$] formally representing 
<I> by matrix multiplication 







\ 


/ 


ai 


\ 


m 


X2 






a2 






\ X„_i 


/ 


V 


an-i 


/ 



(42) 



is lower triangular. 

By combining (|4()|) with the 0(3)-equivariance of $ we obtain the following 
Jacobi vectors 

^ fc-i 

Xfc= Cfc(afc-b(^-i)) = Cfc(afc + ^^^miai) (43) 
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for suitable constants Ck 7^ 0) whose square is determined by the isometry 
condition ii) of namely 

d = =^.l<.-<n-l (44) 
For example, for A; = 1 we have xi = Ci^^i with 



I 1 2 /-2 I 1 2 

|xi| =Ci|ai| 



(m — mi) ■ 



and hence (^f = mim/{ffi — mi). 

Now it is easy to describe all $ satisfying the condition H39|) . namely $ is 
uniquely determined by the choice of sign of the numbers Ck^ ^ — ^ — ''^ ~ ^■ 
There are 2"'"^ such choices, and a natural choice is to take all > 0, namely 
our distinguished $ = <I>o is defined by taking 



C. = y^^^, l<fc<n-l (45) 

Its inverse = will be regarded as our standard Jacobi transformation, 
and hence the standard Jacobi vectors are the columns of the matrix X = 
^'o(X). Then the corresponding lower triangular matrix [^q] = [^o]~^ (see 
(221) )i has entries 

1 i = k 

[*ok = C..<| l<k<i<n-l (46) 

1 <i < k <n-l 

derived from the relation ()43() . On the other hand, the entries of [$o] = (hk) 
are 

1 i = k 

^ l<k<i<n-l (47) 



1< i < k < n-1 



where in each column the entries below the diagonal are identical. Moreover, 
the square sum for each row is given by the simple formula 

lli + - + llk = - l<k<n-l 

m ruk 



Theorem 2.3 The following Jacobi transformation 

^0 = PO ° Pin ° P2n ° ••• ° Pn~l,n ■ Mn ^ M(3, n) 

maps (ai, ...,a„) to (xi, . ..,x„), where x„ = fh~^ ^mjaj is the center of mass 
and the remaining Xj are given by the expressions i4^ ) i4^^ - particular, 
the restriction of to the centered configuration space yields x„ = and it 
coincides with the standard Jacobi transformation 

^'o : Mn M{3,n- 1) : (ai, ...,a„) (xi, . ..,x„_i) 
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The proof follows by an explicit calculation of linear combination of 

ai, ...,afc, which shows that coincides with the vector in ()43() for k < n, and 
furthermore, x„ is the center of mass of the n-configuration X = (ai, ...,a„). 

The above transformation ^'o : Mn M(3, n — 1) turns out, in fact, to be 
identical with the coordinate transformation described in Theorem 1 of Hsiang 
[S], derived by a similar but different "naturality" principle. For convenience, 
explicit formulae (|43j) for the standard Jacobi vectors when n = 3, 4 are listed 
below : 



n = 3 : xi 
n = 4 : xi 

X3 



/ mim 
1712 + rn^ 



ai, X2 



'7712(^2 + 7Ti3^ 



(a2 + 



TTli 



7772 + 7n,3 



ai) 



rrtim 



ai 



r/72 + ?TT'3 + "1-4 
/ 7712(7772 + ?Tl3 + 7774) 

7773 + mi 



(48) 



(a2 + 



777i 



-aij 



' m^{m^ + 777,4) 
7774 



a3 + 



m2 + 7773 + 7774 



(TTliai + 7772a2) 



7773 + 7774 



3 Orthogonal transformation groups on matrix spaces 

By applying a fixed Jacobi transformation ^' : M{3,n — 1), quantities 

and constructions in the centered configuration space M„ are transported to 
its matrix model M(3, tti) with 777 = 77 — 1, as explained in Section 1.4. In this 
setup we may and shall, however, discuss more generally the matrix space 

M = M{d,m) ~ M"^® R'" 

with its Euclidean norm H14|) and orthogonal transformation group 

G = Gi X G2 = 0{d) X 0(777) 

acting by the tensor product Pd® Pmi namely {ip,(p) G G acts on matrices X 
by 

{{ip, X) ^ ipXip ^ (matrix muliplication) (49) 

Physically, we have in mind M as the configuration space for m free bodies (or 
777 + 1 bodies with fixed center of mass) in d-dimensional Euclidean space, and 
the orbit space M = M/Gi is the internal configuration space with the induced 
action of G2 as the internal symmetries. Therefore, for further investigation of 
the topology of M it is natural to analyze its induced G2-orbital structure. Thus 
we shall investigate two successive equivariant systems (Gi,M) and (G2,M), 
where by an equivariant system {K, Y) we mean in general a (compact) Lie 
group K acting on a space Y. We also say {K, y) is a transformation group or 
we simply say y is a X-space, and it is convenient to refer to the orbit space 
construction and the orbit map 

tt:Y ^ Y/K (50) 
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as a K-orhital fibration. In (jSUI) the "fibers" are the orbits, and therefore vr may 
not be the projection of a fiber bundle (or fibration) in the usual sense since in 
general the orbit type is not unique (see below). 

3.1 Compact transformation groups and orbital decomposition 

First we shall recall some basic facts concerning equivariant systems {K, Ai) 
where K is a compact Lie group acting smoothly on a manifold A4, that is, the 
action map K x M ^ M is C°°-smooth. As a good reference on this topic we 
propose, for example, the book Bredonjjj. Then we specialize to the compact 
linear groups, and of particular interest are the natural transformation groups 
on the matrix spaces. 

3.1.1 Basic definitions and constructions 

The space Ai splits into a disjoint union of JT-orbits; these are minimal (or 
homogeneous) i^-spaces in the obvious sense. The isotropy (or stability) groups 
Kp at points p along the same orbit constitute a single conjugacy class (Ka) of 
subgroups representing the type of the orbit. For simplicity we say the orbit is 
of type Ka, and as a -fC-space such an orbit is naturally equivalent to the coset 
space K/Ka- By grouping together orbits of the same type, decomposes as 
a union of orbit strata which are smooth submanifolds Ma, all of which have 
dimension less than Ai except the principal stratum Ai^j, which is an open and 
dense submanifold. The corresponding principal type Ki_j is characterized as 
the unique smallest in the sense that (up to conjugation) C Ka for each a. 
The orbit map 

TT : M ^ M = M/K (51) 

induces a smooth functional structure on the orbit space, in the sense that a 
function f on AA is called smooth if the composition / o vr is a smooth function 
on Ai. Then is a differential space, but it may not be even locally Euclidean 
and hence it is not a manifold in general. However, still it has the nice and rich 
structure of a stratified smooth manifold with the image sets Aia = T^{-M.a) as 
strata, and they are actually smooth manifolds. In this way we may regard Ai 
as a finite union of smooth orbit bundles denoted by 

K/Ka ^ Ma Ma (52) 

whose fibers are those orbits of a fixed type Ka and tTq, is, in fact, the projection 
of an actual fiber bundle. 

For convenience, let us recall the general notion of a (locally trivial) fiber 
bundle, typically illustrated by a sequence like 

F ^ E B (53) 

where vr is the projection map, E (resp. B) is the total (resp. base) space, and F 
denotes the typical fiber, that is, the fibers E}, = 7r^^(6) are homeomorphic to F 
for each b £ B. The simplest example is the product bundle, where E = B x F 
and vr is the obvious projection. An isomorphism between two bundles over B 
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is given by a fiber preserving homeomeorphism (p : Ei ^ E2 between their total 
spaces, and a bundle is trivial if it is isomorphic to a product bundle. The local 
triviality property (which is part of the definition) means that B is covered by 
open sets U such that the portion t[~^{U) over f7 is a trivial bundle, that is, 
7r^^(C/) U >i F . Finally, the bundle is smooth if all the above spaces are 
smooth manifolds and the maps involved are smooth. 

In particular, a JC-orbit bundle like (|52|) is defined by an equivariant system 
{K,E) with a single orbit type, and moreover, B is the orbit space E/K and 
hence the fibers are the i^-orbits. A principal bundle (or principal fibration) is 
the special case 

K ^ E^ E/K = B (54) 

of (|52|) where K acts freely (i.e. all isotropy groups are trivial). 

When K acts by isometries on a Riemannian manifold ^A, there is an in- 
duced Riemannian metric on Aia so that the projection map tTq, in (|52j) is a 
so-called Riemannian submersion. Namely, the tangent map di^a at any point 
p S Ma restricts to an isometry between the horizontal space Tip (i.e. vectors 
perpendicular to the orbit p = K{p)) and the tangent space of Aia at the point 
p. In this way the total orbit space M in 1)51(1 becomes a stratified Riemannian 
space; in particular, the strata are Riemannian manifolds. Moreover, M has 
also a global orbital distance metric which measures the distance between orbits 
in M (c.f. e.g. 0, IHOl)- This metric is certainly determined by its restriction 
to the principal stratum SAu)-, which in turn is derived from the Riemannian 
metric on this stratum. Briefly, the Riemannian manifold determines the 
Riemannian structure on each Ma- 

Next, we turn to the so-called slice theorem which reduces the local study 
of {K^M) to linear representation theory, as follows. The local representation 
of the isotropy group H = Kp at a given point p, denoted by {local)p in 
is the induced linear action (via differentiation) on the tangent space TpM of 
M at p. This splits into the isotropy and slice representation and we write 

TpM =TpK{p) e Hp, {local)p = {Iso)p + {slice)p, (55) 

to indicate that the isotropy representation is the iJ-action on the tangent 
space TpK{p) of the orbit K{p), and the slice representation is the induced 
action on some ff-invariant complementary subspace Tip (which exists since 
H is compact). For example, when M is Riemannian we take Tip to be the 
horizontal space, that is, the orthogonal complement of TpK{p). 

Now, the slice theorem says the orbit K{p) has a i^-invariant tubular neigh- 
borhood which has the structure of a twisted product K Xh Ti-p, consisting of 
equivalence classes [k,v] in X x Tip modulo the relation [k,v] = [kh^^ , hv], for 
k & K,h & H,v £ Tip. The "tube" K TCp has the structure of a vector 
bundle over K{p) ~ K/H with fiber Tip, and moreover, as a -fC-space it has the 
following left action 

k'[k,v] = [k'k,v], k' eK 

and consequently the isotropy groups at points in the tubular neighborhood of 
K{p) looks like 
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In particular, by calculating the slice representation {Kp,Tlp) as the differ- 
ence 

{slice)p = {local)p — {Iso)p (56) 

one can determine the orbit types that occur in the neighborhood. Moreover, 
there is the following isomorphism of orbit spaces 

{K XK, np)/K ~ Hp/Kp (57) 

which enables one to analyse the local smooth structure of orbit spaces by 
successive application of the slice theorem. 

Finally, we recall the notion of a fundamental domain for (K, A^), or equiv- 
alently a cross-section for the orbit map (|51|) . This is a closed subset S C 

which intersects each orbit in a unique point, and consequently there is an 
identification 

S ~ 7W (58) 

In the special case that M is Riemannian and T, is an orthogonal cross-section, 
and hence the submanifold S is also perpendicular to each orbit, the identifica- 
tion (|58|) is in fact an isometry. We remark, however, that fundamental domains 
in the above strict sense exist only in special cases - but fortunately including 
those orthogonal transformation groups that we shall investigate below. 



3.1.2 Compact linear groups on Euclidean spaces 

Now, let us consider pairs (K, V) and the corresponding orbit map 

vr : y ^ F = V/K (59) 

where X is a compact Lie group acting orthogonally on a Euclidean space V. 
Since the orbits are compact it is a well known fact that they can be separated 
by ii'- invariant polynomial functions pi. Namely, we can choose a Hilbert basis 
or "sufficiently many" separating and invariant polynomial functions pi as the 
components of a map 

P = {P,,P2,...,Pn):V 
which induces an embedding 

p:V ^ p{V) C (60) 

and hence identifies the orbit space with its image p(T^). The latter is a semi- 
algebraic subset, that is, it is defined by polynomial identities and inequalities 
in M^. 

Concerning smoothness, we recall that V has, on the one hand, the structure 
of a differential space such that a function / on F is " smooth" if and only if the 
composition /ovr is a smooth function on V. On the other hand, there is another 
approach defining a function on p{V) to be "smooth" if it is the restriction 
of a smooth function on M^. Fortunately, the two notions of smoothness are 
identical, and therefore we may regard the bijective correspondence p in (|6()j) to 
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be a diffeomorphism, cf. Schwarz . Consequently, / is " smooth" if and only 
if/o7r = /op for some smooth function / on R-^. Note, however, the orbital 
distance metric on y ~ p(^) is certainly not induced from the Euclidean space 
in PH). 

In the Euclidean space V size and shape are quantities preserved by or- 
thogonal transformations. The norm p = |X| is a size function which measures 
the distance from the origin O, namely the arc- length along rays emanating 
from O. A nonzero vector can be uniquely scaled to unit size by a homothety 
transformation X — > cX (c > 0), and by expressing the Euclidean metric ds'^ 
of V in polar coordinates we actually describe V as a Riemannian cone over its 
unit sphere : 

{V\de^) C {V,ds'^), ds'^ = dp^ + p'^de'^ (61) 

where dO'^ denotes the induced spherical metric on V^. 

On the other hand, the JT-orbit space V and its subspace V* = V^/K are 
stratified Riemannian spaces as explained above. Moreover, since the norm 
function p is i^T-invariant it is also a size function on V and then V* = [p = 1) 
becomes the unit distance "sphere". Furtermore, in analogy with (|6H) V inherits 
the structure of a stratified Riemannian cone over its "sphere" V* : 

{V*,da^) C {V,df), ds'^ = dp^ + p^da^ (62) 

where da"^ denotes the orbital distance metric induced from dO"^. In this cone 
the rays emanating from the cone vertex (or base point), still denoted by O, are 
also geodesies for the metric ds^, and homothety transformations move points 
7^ O along these rays. 



3.1.3 The configuration space for m-body systems in d-space 

Henceforth, we take V to be the matrix space M = M(d, m) with the orthogonal 
action of G = 0(d) x 0{m) as in ()49() . By regarding M as the configuration 
space for an m-body system in W^, in analogy with the "physical world" case 
d = 3 described in Section 1, Gi = 0{d) and G2 = 0{m) play the roles 
of congruence and internal symmetry groups respectively. Then there is the 
following two-step orbital decomposition of stratified Riemannian spaces 

M M M 

-ar err a 

where the internal configuration space M = M(d, m) and its symmetry group 
constitute a natural equivariant system {G2^M). Again, the metric ()62() on M 
will be referred to as the kinematic metric^ since in the physical case d < 3 it 
is, in fact, representing the internal kinetic energy 

f = -df = T -T"^ = -ds^ - T'^, 
2 2 

namely the total kinetic energy T minus the purely rotational kinetic energy 
T'^ of the n-body system, see e.g. [30] . 
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The internal symmetry orbit space M/G2, which is a cone in M^, /i = 
min{d, m}, is actually a Weyl chamber of Cartan type B^. Indeed, the G- 
representation on M is the isotropy representation of the symmetric space 0{d+ 
m)/0{d) X 0(m), which is of type B^, and as such it is well understood in 
terms of the classical invariant theory of matrices and canonical forms modulo 
orthogonal transformation groups. 

Let be the unit sphere of M. In view of the cone structure (|61|1 . (|62|) . 
we may faithfully replace by its unit distance "sphere" version 

..I M* Ml 

^ M* = — > — = — = M** (64) 

where M* = M* (d, m) will be referred to as the shape space. Its G2-orbit space 
M** is a spherical simplex and it is explicitly described in (|85j) below. 

The orbit types of the above transformation groups will be described in 
detail. In principle, the topological structure of the shape space M* can be 
reconstructed and analyzed as the union of well understood orbit bundles, but 
in this paper we shall only describe the global topology of M* in those few 
cases where it is a msnifold. In particular, the cases d = m — \ or d = m are 
rather simple since (G2, M*) is essentially isomorphic to a linear model, that is, 
a sphere or disk D^~^^ with an orthogonal action of G2, see Section 3.4.3. 
Before turning to these topics, let us have a closer look at some important 
classical matrix spaces and their natural transformation groups. 



3.2 Transformations of the space of symmetric matrices 

For d < m, consider the following sets of m x m-matrices 

Sym'^{m)<(i C Sym^{m) C Sym{m) C M{m,m) 

U ~ U U (65) 

B+{d) ^ B+{m)<d C B+{m) C B{m) C Diag{m) ~ M™ 

where in the first row Sym{m) is the space of symmetric matrices, Sym~^{m) 
is the convex cone of positive semidefinite matrices, and Sym~^{m)<(i is the 
subcone of matrices of rank at most d. In the second row Diag{m) is the space 
of diagonal matrices diag(Xi, Am), containing the cone 

B{m) = {Y = diag{Xi,X2, .■.,Xm); Ai > A2 > . .. > Am} , (66) 

and the subcones B'^{m) and i?^(?7i)<^ are given by the additional constraint 
Am ^ 0, and Aj = for all i > d, respectively. 
With the usual inner product on M{m, m) 

Yi-Y2 = trace{YlY2), cf. (dU 

where y* means the transpose of Y, 0{m) acts orthogonally by conjugation 

{if, Y) ipYif-^, if E 0(m), Y G M(m, m) (67) 

and the subsets in the first row of (jHSJ) are clearly invariant. Thus the transfor- 
mation group {0{m), Sym{m)) is the second symmetric tensor product S'^p^ 
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of the standard representation = (0(m),M™'), and let us recall the classical 
result about diagonalization of symmetric matrices. Namely, every symmetric 
matrix is conjugate to a unique matrix - its canonical form - in the cone B{m), 
which is therefore a fundamental domain (or cross section) of the orbit map 

^ , , Sym{m) M™ . ^ ^ . ^ 
Sym{m) ——. — r— ~ — ~ B(m) C Sym(m) (68) 
U(rn) Dm 

As indicated, B(m) is also a fundamental domain (or cross section) for {Sm, K"*), 
where 5m is the symmetric group acting on M"^ by permuting the coordinates 
A.. 

Fundamental domains for the 0(m)-action restricted to invariant subspaces 
of Sym{m) are the corresponding subsets of B{m) in the second row of (|65() . 
consequently 

As indicated, we recognize B~^{k) as a Weyl chamber for the compact connected 
Lie group SO{2k + 1), namely a fundamental domain of the canonical repre- 
sentation of the Weyl group Bj^ D Sk acting on M.^ as a group generated by 
reflections. 

Note that Sym[m) is not an irreducible 0(m)-space. In fact, there is the 
orthogonal and invariant decomposition 

Sym{m) = Sym^{m) ® m} (70) 

where the trivial summand is spanned by the identity matrix, and the first 
summand, which consists of matrices of trace zero, gives the irreducible rep- 
resentation S'^Pm — 1 with fundamental domain B^{m) C B{m) defined by 
EA. = 0: 

Here we have also recognized B^im) as a fundamental domain of the canonical 
representation of the Weyl group Am-i = Sm of SU{m). 

In addition to its Euclidean norm where \Y\ = y^trace{Y'^), the positive 
cone Sym'^{m) has another 0(m)-invariant size function p*(Y) = trace{Y), 
with corresponding "unit sphere" 

Sym*{m) = {y G Sym^{m) ; trace{Y) = l} (72) 

For later reference, we note that its intersection with the Euclidean unit sphere 
in Sym{m) is the subset 

Sym*{rn) n Sym{m)^ = Sym* (711)1 = {Y £ Sym*(m) ; rk{Y) = 1} (73) 

Let Y G B{m). The tangent space at Y of the 0(m)-orbit through Y in 
Sym{m) is the subspace 

[so(m), Y] = {SY — YS ; S is skew-symmetric} C Sym{m), 
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and it is perpendicular to B{m) since the inner product with any Y € B(m) is 

{SY - YS) ■ Y' = trace{{SY - YS)Y') = 
A direct consequence of this observation is as follows : 

Remark 3.1 The fundamental domain B{m) of {0{m), Sym{m)) is an or- 
thogonal cross section, in the sense that it is perpendicular to the O{m)-orbits. 
Therefore, the orbit space of Sym[m) (with the orbital distance metric) is iso- 
metric to B{m). Consequently, orbit spaces of various 0{m)- invariant sub- 
spaces of Sym{m) in i^,5)) are isometric to the corresponding subsets of B{m). 

Finally, let us determine the orbit types of the 0(m)-action on symmetric 
matrices. To this end, consider a matrix Y G B{m) and let the "strings" of 
equal entries Aj in have length mi,m2, so that 

m = mi + m2 + ... + mp {rrii > 0, p > 0) 

is a partition of m and p is the number of different X'^s. It is easy to see that 
the isotropy group at Y is the "block" orthogonal matrix group 

0{mi,m2, mp) = 0{mi) x 0{m2) x ... x 0{mp) C 0(m) (74) 

and therefore the partition function 7r(m) of m enumerates the different orbit 
types. The orbits are in fact connected and therefore they are also SO{m)- 
orbits, that is, 

Sym(m) Sym(m) ^, . 

— - = — - ~ H( m) 

SO{m) 0(m) " ^ ' 

3.3 Algebraic realization of orbit spaces and orbital stratifica- 
tion 

With the results from the previous subsection as a basis, the construction of the 
orbit spaces in (|63j) or (|64|) as semialgebraic subsets is based upon the following 
two simple but fundamental properties of matrices : 

• The polynomial map 

p:M{d,m) Sym{m), X^Y = X*X (75) 

preserves the matrix rank, and its image consists of all positive semidef- 
inite symmetric matrices Y of rank < min{(i, m}. Moreover, for ij) G 
0{d),^e 0{m), 

piijXip-') = ^p{X)^-' =(^y^-i (76) 

and the entries of Y are the inner products Xj • xj of the column vectors 
Xj of X. 

• Consider a finite collection of vectors Xj in a Euclidean space W^. Then 
the collection is uniquely determined, up to 0(d)-congruence, by the inner 
products Xj • Xj. For a proof we refer to Weyl|3X1, page 52. 
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The map p in ((TSJ) is clearly constant on each 0((i)-orbit in M = M{d, m) 
and induces a surjective map 

p : M{d,m) Sym~^{m)^d C Sym{m) (77) 

which is, moreover, injective by the above mentioned result of Weyl. In other 
words, the components of p, namely the collection of inner products Xj • Xj, 
constitute a Hilbert basis in the sense of Section 3.1.1. Consequently, p is a 
diffeomorphism between differential spaces which identifies M = M/0{m) with 
the set of positive semidefinite symmetric m x m-matrices of rank at most d. 
In particular, 

M{d,m) = M{m,m) ^ Syrn'^{m) if d>m, (78) 

so we will henceforth assume without loss of generality that d < m. Thus we 
have explained why all orbit spaces involved are realizable as appropriate spaces 
of symmetric matrices, amenable to the setting in Section 3.2. 



3.3.1 Cross sections and canonical foms 

In general, the inclusions M{k,m) C M{k + l,m), defined by taking the last 
row to be zero, induce the increasing filtration 

M(l,m) C M(2,m) C ... C M(m,m) ~ S'ym+(m) (79) 

and this coincides with the matrix rank filtration of the positive cone Sym^^m), 
namely for k + i < m 

M{k, m) = M{k + i, m)<k ~ Sym'^ {m)<k 

On the other hand, G2 = 0{m) acts on the spaces in (|77() . The action 
on M = M(d,m) is induced from the action ()49() . and G2 acts on Sym{m) 
by conjugation 1)671) . From the equivariance property (|76|) it follows that these 
actions commute with the map p, that is, p is an isomorphism of G2-spaces and 
hence induces a diffeomorphism between G2-orbit spaces 

P •7T = 7^ n - C Sym+{m)<d - M (80) 

The cross section B^[d) C M in ()8U|). being transversal to the G2-orbits in M, 
further lifts to a cross section M~^{d) of the composite orbit map M M/G, 
as follows : 

B^d) M^d) C M : I ^ = ^^«^(n,r2, ... r,) 

where X = {xij) has xa = ri and zero entries otherwise. Thus we have also 
established the diagonalization procedure saying that every matrix in M = 
M(d, m) can be transformed by the action (|49)) to a unique matrix - its canonical 
form - in the subset M^{d). As a (simplicial) Euclidean cone in the {rj}- 
coordinate space W^, M^{d) is also the Weyl chamber of the Weyl group 5^, 
cf. (|ni). 
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3.3.2 Orthogonality of cross sections 

In view of Remark 13.11 the cone (d) in l|81jl is orthogonal to the G-orbits 
and is therefore an orthogonal cross section for (G, M). Hence, we have an 
isometry M/G ~ M~^{d) of Euclidean cones with the metric 

i=l i=l * 

Here {rj} from (|81() and {Aj} from H66() are coordinate systems for the cone and 
are related by the polynomial map p in (jTSjl , which restricts to a diffeomorphism 

p: M+{d) ^ B+{d) C Sym+{m), n ^ rf = \i (83) 

Thus, the G-orbit of a given matrix X G M = M{d, m) has coordinates Aj 
interpreted as the eigenvalues of X^X. In the special case that d = m and X 
is symmetric, the numbers are the absolute values of the eigenvalues of X. 

Remark 3.2 At this point, observe that the internal configuration space M = 
M{d,m) has two metrics, namely the kinematic metric (i.e. the 0{d)-orbital 
distance metric) and the induced "fake" metric as a subset Sym~^{m)<^d of the 
Euclidean space Sym(m), cf. i7T\ ). For both metrics the symmetry group 0{m) 
is actually an isometric transformation group with B^[d) as a fundamental 
domain (i.e. cross section). In the "fake" metric B~^{d) inherits the Euclidean 
metric dA? and is by Remark \S.l\ an orthogonal cross section, whereas in the 
kinematic (but still Euclidean) metric i8^) of M the orthogonality property of 
B^{d) fails. See also Remark \cl.4\ below. 

3.3.3 Rank and subrank stratification of matrix spaces 

The subrank stratification of the matrix space M = M{d, m) is a natural refine- 
ment of the usual rank stratification. The latter coincides with the orbit type 
stratification of the action of Gi(or G2) on M and the mentioned refinement 
coincides with the orbit type stratification of the action of the full group 

G = Gi X G2 = 0{d) X 0(m) 

In order to describe the combinatorial structure involved we turn to the two- 
step orbital decomposition (jS^, where the first orbit space is the shape space 

M* = M*{d,m) ~ Sym*{m)<d , cf. JZl (84) 
and the final G-orbit space is 

M** = M**{d,m) = M*/G2 = M^/G 

which we may identify with the following fundamental domain in the sphere 

A'^-i = M+{d) n : n > rs > ... > > 0, ^ = 1 (85) 
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This is a spherical simplex whose structure, indeed, identifies it with the spher- 
ical Weyl chamber of type Bd (see (|69() ): in particular, it is homeomorphic to 
the closed disk D'^~^. 

Thus we also conclude that the above simplex H85|) is a cross section for the 
transformation group (G, M^), and hence any matrix in is G-equivalent to 
a unique matrix - its canonical form - in the simplex. On the other hand, a 
typical matrix of rank k in (|85j) may be characterized by the following numerical 
data 



^ { 'ri=rk^> rfc^+i = rfc^+fca > ^fci+fc2+i = •••• = '^fc > = r^+i , , 
I 1 < ELi ki = k<d, ki>0, p>0, ^"""^ 

and then its subrank is defined to be the corresponding (unordered) partition 
of k into p positive integers 

K = {ki,k2, ...jkp) , |k| = ki+ ... + kp = k (87) 

which records the length of the strings of equal numbers in (|86|1 . 

We label the rank (resp. subrank) stratum by the subscript k (resp. k) of 
the corresponding sets. Thus X belongs to if and only if its Gi-orbit has 
the type 0{d — k) (cf. also (|inn|) below). Moreover, by (|85|) the image of 
in M** is the "semi-open" simplicial disk 

Mr ^ A^-i C : ri > r2 > ... > > r^+i = 0, J^r^ = 1, 

which is subdivided into its various subrank strata 

Ml* ~ A^~^ C A^-^ C A'^-i ~ M** (88) 

In particular, the principal stratum has subrank kq = (1, 1, 1) with |ko| = d, 
and M** ~ A^~^ is the interior of the spherical simplex A'^"^ and is therefore 
homeomorphic to M^~^. 

To see why the subrank stratification actually coincides with the G-orbit 
type stratification, consider the isotropy group Gx of the matrix in (|5H|) . using 
the following notation for groups : 

A.H = diag{H x H) ~ H (diagonal embedding) 
0{ki, .., kp) = 0{ki) X ... X 0{kp) c 0{k) (89) 
G{ki, .., kp) = 0{d -k)x AO{ki, .., kp) x 0{m - k) c G 

Then it is not difficult to show by direct calculations with matrices that 

Gx = {(91,92) eO{d) X 0{m) ; g^Xg^^ = X} = G{k,, ..,kp) (90) 

(cf. Table 1, #10, in [17j), and hence the isotropy types uniquely characterize 
the subrank strata, as claimed. 

We are particularly interested in the stratification of M*, but now there 
are two natural options, namely the induced subrank stratification and the G2- 
orbit type stratification. Our next claim, however, is that they are identical. 
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and hence one can describe M* is a union of G2-orbit bundles lying over the 
various strata (|55)) . 

To calculate the G2-orbit types and how they correspond to the subrank 
strata we proceed as follows. Let X* £ M* be the image of a matrix X of 
rank k and consider the "large" orbit G{X) C M. As a Gi-space all the orbits 
of G{X) have the same type as the "small" orbit Gi{X). Moreover, the orbit 
space G{X)/Gi C M* coincides with the G2-orbit through X* , namely the set 
G2{X*). Thus, the "large" G-orbit through X has the structure of a Gi-orbit 
bundle, and in terms of homogeneous spaces (i.e. coset spaces of groups) this 
fiber bundle can be described as follows : 

Gi{X) ^ G{X) ^ G2{X*) 

M II II (91) 

0<,d) G . 0{m) 

0{d-k) G(ki,k2, -^kp) 0{ki,k2, ..,kp)xO(m-k) 

The base space G2{X*) of the fiber bundle is in turn a fiber of another bundle, 
namely the multistratum M* which as a G2-orbit bundle fibers as follows 

'^^"'^ M* M** ~ A^-i (92) 



0(^1, k2, .., kp) X 0{m — k) 
In this way the subrank stratification of the shape space 

d d 
M*{d,m) = \jMl=\J U M: 

k=l A:=l|K|=fc 

has strata which are fiber bundles of type (|92() . and there are althogether 

7r(l)+7r(2) + ... + 7r(d) (93) 

such bundles (or strata), where 7r(A;) is the partition function of k. The principal 
stratum M*^ constitutes an open and dense subset of M*, and the principal 
orbit bundle 

'M*-M:>M'^-i (94) 



O(l)'^ X 0{m-d) 

is actually trivial (e.g. since the base space is contractible) . The orbit type, 
either for the G2-action on M* or the G-action on M^, for a stratum of subrank 
K = {ki, k2, kp) can be read off from (|91() . 

Example 3.3 d = 3 : M(3,n— 1) ~ M„ is the centered configuration space for 
n-body systems in 3-space. See Section 3.4-4 for a brief description of the shape 
spaces M* when n is small. In general, for n > 3 their G2- orbit space 

M* 

^ = Ar*(3,n-1) ~ 
G2 

is the spherical triangle on the sphere + + r| = 1 with vertices 

^ = (1,0,0), 5 = (1/^2,1/^/2,0), C = (1/^3, 1/V3, 1/^) 
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On the circle = lies the closed circular arc of length 7r/4 

[A, B]:ri>r2>0 

which consists of those points of rank < 2. It has three strata, namely the vertex 
A (of rank 1), the vertex B and the open arc {A,B) (both of rank 2). However, 
the open arcs 

{A, C) : ri > r2 = r-s > 0, {B, C) : n = r2 > rs > 

constitute the same stratum of subrank k = (2,1). Thus, although the triangle 
decomposes into 7 = 3 + 3 + 1 strata components, namely vertices, open edges 
and the interior of A^, there are only 6 multistrata (since one of them has 2 
components), in accordance with the enumeration formula \9S\) . 

Remark 3.4 In general, the kinematic metric on the shape space M* is uniquely 
determined by its restriction to the principal stratum, which is a Riemannian 
manifold (M*^,da'^). There is a general procedure, in the setting of equivariant 
differential geometry, for the calculation (or description) of the Riemannian 
connection on such a principal orbit bundle, see for example The simplest 
case d = 3, as in the above example, is analyzed in HsiangJ^, but similar cal- 
culations can also be done for d > 3. We shall, however, leave the geometric 
issues here and concentrate on the topological structures in the sequel. 

3.4 Topology of the shape space and related spaces 

For fixed m the largest shape space (with d = m) is 

M*{m,m) ~ Sym*{m) (95) 

and by (|H1|) it contains any other shape space 

M*{d,m) = M* {m,m)<_(i ~ Sym* {m)<_d 

as the union of those strata of matrix rank at most d. Therefore, we shah refer 
to the space (|95|) as the m-universal shape space. Thus there is the increasing 
and 0(m)-invariant rank filtration 

M*{l,m) C M*(2,m) C ... C M*{m,m) (96) 

and corresponding 0(m)-orbit spaces ~ A*^^^, 1 < d < m, 

M**(l,m) C M**{2,m) C ... C M**{m,m) = A"""^ (97) 

where the first three spaces in the chain (|97j) are, indeed, the inclusions of 
simplices 

{A} = A° C [A, B] = A^C [A, B, C] = A^ 
from Example 13.31 
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By inspection the chain in ()96() starts with 

M*{l,m) = 5™-VO(l) = MP™-^ 

M*(2,m) = 52™-VO(2) = CP"'-^/Z2 (98) 

where 

are the real and complex the projective (m-l)-space, and Z2 = 0{2) / S0{2) 
acts on CP™""^ by complex conjugation with MP™""^ as fixed point set. In 
particular, the ver tex A° = M**(l, m) of the simplex A^-^ in (j!?!) is the single 
0(m)-orbit M*(l,m), with the topology of a real projective space as noted 
above. 

The topology of the next shape space M*(2, m) is more difficult to describe. 
One approach is to utilize the fact that it is a 0(m)-space of cohomogeneity 
one, that is, with one-dimensional orbit space. This holds since 

^^^^ = M**{2,m)^A' = [A,B] 
(J[m) 

is an interval. We also know the orbit types (i^i), {K2), (H) corresponding 
to the three strata A, B and {A, B) respectively. Therefore, the space can be 
described topologically in terms of its orbit types by the construction (cf. e.g. 
1291, Vol.1, Chap. IV) known as the equivariant union 

M*{2, m) = M(7ri) U M{tt2) 

of the mapping sylinders M(7rj) of the canonical projections vrj : 0{m)/H — > 
0{m)/Ki, where in our case 

Ki = 0(1) X 0{m - 1), K2 = 0(2) X 0{m - 2), (99) 
H = KinK2 = 0{lf X 0(m - 2) 

The space M*(2,2) is a 2-disk and M*(2,3) is actually homeomorphic to 
S^, see p()4|) and Section 3.4.4). For m> 3 M*{2,m) fails to be a manifold in 
a neighborhood of the singular 0(m)-orbit 

A = 0{m)/Ki = M*{l,m) ~ MP'"-\ 

but two copies of M*{2,m) glued together along A yields, indeed, the dif- 
ferentiable manifold S"^^'^ /S0{2) = CP™~^. Below we wih return to this 
construction and further investigate the topology of the spaces M* (d, m) for 
2 < d<m. 

3.4.1 Local and global topology 

By referring to Section 3.1 concerning compact transformation groups and the 
slice theorem, let us first recall a nice property of regular representations of 
classical groups such as 0{d). Namely, we consider the linear action of 0{d) 
on by some representation of type $ = mp^ + Tg, N = md -\- q, which is m 
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copies of the standard representation plus a g-dimensional trivial summand 
Tq. Then the orbit types constitute the following "connected" string of regular 
subgroups 

0{k) ; d-m <k<d (100) 

where 0{k) = 1 if A; < 0. To calculate the slice representation <I>fc at a point 
with isotropy group 0{k) we calculate the difference as in ()56p . with the local 
representation equal to ^ restricted to 0{k), that is, mpj. + {trivial) . Moreover, 
the isotropy representation Iso{d,k) of 0{k) is the "linearized" action on the 
Stiefel manifold 0(d) /0{k) = SO{d) / SO{k) at the base point, and a simple cal- 
culation (involving the adjoint representation of 0{d)) yields the representation 
{d — k)pj. + (trivial), consequently 

$fc = ^\o{k) ~ Iso(d, k) = (m — d + k)pf. (mod trivial ) (101) 

From this we conclude that each slice representation inherits the regularity 
property of the original action 

Now we turn to the "left side" action of 0(d) on the matrix space M = 
M(d, m) ~ M."^'^ , which is just the regular representation ^ = mp^. However, 
we shall rather consider the restricted action on the unit sphere M^, noting 
that the only change in the above calculations is that the local representations 
lose a trivial summand ri. The orbit types still constitute the string (jlOOl) and 
also holds, except that k < d \n since 0(d) has no fixed point on 

the sphere M^. 

Lemma 3.5 (i) The m-universal shape space M*(m,m) ~ Sym*(m) is a man- 
ifold with boundary 

dM*{m,m) = M*(?tt- — 1, m) ~ «S'ym*(m)<m-i 

(ii) For p > 0, M*(d,d + p) is a manifold (and with no boundary) if and 
only if p = 1 or d = 1. 

Proof. By ()57() . the local topology of M*(m,m) around an 0(m)-orbit of 
type 0(k) is the topology of the orbit space of the slice representation (0(k), ^k) 
in M^, with m = d in Hence, part (i) of the lemma follows by induction 

on m > 1, since at the initial step m = 1 the orbit space of (0(l),pi) is the 
half-hne M.X = [0, oo) . 

For the proof of (ii) we may assume d > 1, and the slice representation ()101() 
of 0(k) is = (k + p)Pk (mod trivial). By induction the proof reduces to the 
crucial singular case (0(1),(1 -\-p)pi), where 0(1) acts by inversion v — > —v, 
and here the orbit space is the cone = C(M.PP) over the real projective 
p-dimensional space. This is a manifold if and only li p = 1, in which case the 
orbit space is homeomorphic to the Euclidean plane 

^ = C(RP^) « C(S^) = (102) 
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Next, let us actually determine the topological type of the compact manifold 
M*{m,m) ~ Sym*{m) and its boundary. By the 0(m)-action on M{m,m) 
every X = (xij) can be mapped to the subset 

T^(m) C M{m,m) 

of upper triangular matrices with diagonal entries xu > 0, and therefore the 
restriction of the polynomial map in ()75() 

p : T+(m) ^ 5ym+(m) (103) 

is still surjective. The subset T^{m)m of matrices X with all xa > 0, i.e. of 
maximal rank m, is clearly diffeomorphic to the Euclidean space M™(™-+^)/2^ 
and it is easy to verify that different matrices X lie on different 0(m)-orbits. 
Therefore the following lemma must hold. 

Lemma 3.6 The polynomial map in \1US\) further restricts to a diffeomor- 
phism 

p : r+(m)„ ~ Sym+{mlrn ^ M"^(™+i)/2 

with the set of positive definite symmetric matrices. Hence, the interior of 
Sym*{m) is the open "sphere- octant" 

Sym*{m)m ^ T+{m)m nM^{m,m) : = l,xu>0 

which is diffeomorphic to K"t-('"+i)/2-i _ 

We know from the previous two lemmas that Sym*{m) is a compact man- 
ifold with boundary and its interior is an open disk. Using some manifold 
theory this information actually suffices to conclude that Sym*{m) is a closed 
disk. For example, one may apply the so-called collar neighborhood theorem 
and the h-cobordism theorem, cf. Milnor[^). Anyhow, we have the topological 
types 

M*(m, m) ~ Sym*{m) w D™('»+i)/2-i (104) 
M*(m - l,m) ~ Sym*{m)<m-i « 5'"(™+i)/2-2 

and taking the cone over these spaces yields the topological types 
M(m, m) « M^^^+i^/^-i ^ [q, 00) (half-space) 
M{m - 1, m) ~ dM{m, m) « m™(™+i)/2-i (Euclidean space) 

Remark 3.7 The linear model construction in Section 3.4-3 provides another 
proof of the homeomorphisms \1U4\I - 
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3.4.2 Branched coverings and generalized Hopf fibrations 



The usage of SO{d) rather than 0{d) as the congruence group acting on M{d, m) 
may lead to a different " shape space" , and we shah explain this distinction be- 
low. First, observe that for k > 1 we have connected and hence equal orbits 

0{d)/0{k) ~ SO{d)/SO{k) 

and consequently the orbit spaces of SO{d) and 0{d) will be identical as long 
as d > m, see (|100|) . However, due to (|78]) we have been assuming d < m and 
therefore A; = also occurs in the string (|1UU() . so the two orbit spaces cannot 
be identical. In order to explain the difference the following discussion will be 
helpful. 

Consider a space Q with a given involution a (i.e. a transformation of order 
two) and hence we shall regard Z2 = {Id, a} as a transformation group on Q. 
We denote the orbit space of Z2 by Q and the fixed point set by S. We also 
assume E 7^ and regard it as a subset of both Q and Q. Then the orbit map 

n:Q^^=Q (105) 

is an example of a ramified double covering which is ramified along S. Con- 
versely, starting from Q and T, we may reconstruct Q as the double of Q 

2Q = QU^Qc^Q 

by taking two copies of Q and identify (or glue together) their "singular" set 
E. However, the pair {Q, S) does not always lead to a unique double space, so 
we shall rather have in mind a specified ramified covering as in H1U5|) . Thus, we 
define for d < m 

2M*{d, m) = M*{d, m) UM-(d-i,m) M*{d, m) (106) 
and refer to the orbit map 

2M*(<i,m):.-^ -^-^ = M^d,m) (107) 

of the action of 0{d) / SO{d) = Z2, which is, indeed, a double covering ramified 
along the fixed point set M*{d — 1, m) of Z2. 

Lemma 3.8 For p > and d> 1, the space 2M*{d,d + p) is a manifold (and 
with no boundary) if and only if p = or d = 2. 

Proof. For d = 2 the 50(2)-orbit space is the complex projective space 
CP^"*"^, so let us assume d > 2. Using again the slice theorem and an inductive 
argument, the crucial case will be the topology of the orbit space of {S0{2), {p-\- 
2)^2)) namely the cone 

^2{2+p) c2(2+p)-l 

50(2) = ^(^0(2r> = '^'^^ ' 
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over the complex projective space CP^^^. It is well known that this is a man- 
ifold (i.e. locally Euclidean) if and only if p = 0, in which case it is a cone 
homeomorphic to Euclidean 3-space 

C(CP^) w ^(5^) = M.^ (108) 



Lemma 3.9 The induced smooth functional structure on the cone \10^) (resp. 
mUSp ) is a refinement of the standard differentiable structure o/M^ (resp. W^). 
The structures are identical away from the origin (cone vertex), where the cone 
fails to be a smooth manifold. 

Proof. Consider the case ()1U8() . namely the orbit space of (50(2), 2/?2) 
acting on x M^. To describe it algebraically, let u, v denote vectors in 
and define the polynomial map 

p = {X, Y, Z, VF) : IR2 X m2 ^ , see ^ 

where X = juj'^ — |v|^ ,Y = 2u • v, Z = 2(u x v) • k, W = |u|^ + |v|^, which 
identifies the orbit space with the cone in 4-space 

x'^ + y'^ + z"^ = , W>0 

By projecting the cone onto the coordinate {X, Y, Z)-space we identify it with 
as in p()8|) . with (standard) coordinate functions X, Y, Z which are certainly 
smooth in the induced orbital functional structure. However, at the origin the 
function / = \/X'^ + Y'^ + Z'^ is not smooth with respect to the standard struc- 
ture, although its composition with the orbit map is by our construction the 
smooth function W. The case (|1()2|) is similar but simpler, using p = {X, Y, W) 
and {u,v) e x ■ 

According to (|104|) and (|107j) we have constructions of orbital fibrations of 
spheres over spheres involving two types of orthogonal transformation groups, 
as follows : 

Case 1 : Consider the regular representation {0{m — l),mp^_i), acting on 
the unit sphere M^{m — l,m) = ^l*^-!)™-! with orbit space as in the 
second line of (|1U4I) . 

0(01—1)771—1 

^(m-l)m-l ^ ^ ^ ^m(m+l)/2-2^ m = 2,3,4,... (109) 

0[m — 1) 

Here the initial case m = 2 is the trivial Hopf fibration or double covering 

c-l 

0(1) ^ RP^ = ~ (110) 

Case 2 : Consider the regular representation {SO{m),mp^), acting on the 
unit sphere M^(m,m) = ~^ . The orbit space coincides with the 
"double" construction (|10tj|) applied to the first line of (|104p 

= 2M*(m,m) ^ 20^^"^+'^^-' = S^^^+^^'^-^ 

bO(m) 
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and there is the orbit map 



^m(m+l)/2-l^ m = 2,3,4,.... (Ill) 



SO{m] 

which for m = 2 is the well known Hopf fibration 

SO{2) ^ 53 CP' = (112) 

Remark 3.10 With the induced differential structure the above quotient spheres 
are, of course, stratified differentiahle manifolds ( cf. Section 3.1) which are also 
locally Euclidean spaces since they are spheres in the topological sense. In fact, 
for m = 2 they are the standard sphere (i.e. S' or S'^) in the differentiahle 
sense as well. However, for m> 2 there are more than one orbit type, and with 
the induced differential structure they are actually not differentiable manifolds. 
However, this artificial situation can be remedied by "relaxing" the differential 
structure so that the quotient sphere becomes the standard sphere while the map 
IT is still differentiable, see Section 3.4.3. 

Remark 3.11 The above orbital fibrations are natural maps between spheres 
with specific properties, with possible future applications in physics. Indeed, the 
Hopf fibration Mli^) is a simple special case with several well known applica- 
tions. For example, it describes the geometry of a magnetic monopole, and 
Dirac made the major discovery that the fibration could explain (in a modern 
language) the quantization of electric charge. Here is the unit sphere of 
= C^, the action of 5*0(2) ~ U{1) extends to by complex scalar multi- 
plication, and the point magnetic source is at the origin of C? /U{\) ~ M^. The 
4-space in this setting is also referred to as the Kaluza-Klein model of the 
Dirac monopole. 

Arnold pi has considered the three special cases m = 2, 3, 5 of (|lfl9j) from a 
different viewpoint. Namely, they fit into the following unifying pattern with 
projective spaces over R, C or EI as intermediate quotient spaces : 

i) S' MP^ ~ S' 
ii) ^ CPV0(1) « (113) 

Hi) ^ M ^ M P^/ (Anffl, Z2) ^ S^'^ 

where 0(1) — acts by (complex) conjugation and AutM. ~ 5*0(3) is the 
automorphism group of EI (quaternions). The last two cases are two-step orbital 
fibrations defined by groups H ^ K ^ K/H^ that is, K acts on the sphere 
on the left side in (|113j) . and the intermediate (projective) space is the orbit 
space of the normal subgroup H, with the induced action of K/H. The groups 
corresponding to the last two cases of (|113|) are as follows 

ii) 30(2) 0(2) ^ 0(1), Hi) Sp{l) 0(4) ^ 0(3) 
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where 0(3) = SO{3) x Z2 — (^utH, Z2). However, from this viewpoint it is 
also natural to include the cases m = 2, 4 of where there is no group Z2. 

Case m = 2 is the Hopf fibration H112jl and the new case m = 4 reads 



with the corresponding groups -5^(1) 50(4) SO{3). 
3.4.3 The m-universal linear model 

We shall describe another prperty of the m-universal shape space M*(m,m) ~ 
Sym*{m), which by (|104() is known to be a disk whose boundary sphere is 
M*{m — l,m). In fact, what is remarkable is that the equivariant system 
{0{m), Sym* (m)) has a linear model, namely it "resembles" closely a Euclidean 
disk with 0{m) acting as an orthogonal transformation group. To make this 
correspondence precise, let us first inquire what is the appropriate linear model. 

Recall from Section 3.2 that 0{m) acts orthogonally on the space Sym{m) 
of symmetric matrices of dimension m, namely by the symmetric tensor product 
representation S'^p^^. Let SymP{m) be the subspace of matrices of trace zero, 
where 0{m) acts by the irreducible representation S'^p^ — 1, and let 

D^{m) C Sym^{m) cf. ^ 

be the unit disk centered at the origin. Then we make the following definition 
which will be justified below. 

Definition 3.12 The above equivariant system {0{rn),D^{m)) is the linear 
model of {0{m), M*{m, m)), and the restriction to the boundary sphere dD^{m) 
is the linear model of dM* {m,m) = M*{m — l,m). 

Now, we shall construct a 1-1 correspondence between Sym*{m) and the 
disk D^im) which is smooth and 0{m)- equivariant, that is, the map commutes 
with the action of 0{m). To each matrix Y G Sym*{m), namely a positive 
semidefinite matrix with trace 1, we associate the matrix 



H p3 ^ M P^/Autm ^ S' 



(114) 



Yq = — (Y Id) G Sym (m), where r. 




and note that Yq is perpendicular to the identity Id and 




4-(irace(y2) - — ) < ■^(trace(Y)'^ - — ) = 
r4 m m 



Therefore, the 0(m)-equivariant afhne transformation 



Sym{m) Sym{m) : Y — (Y Id) = Yq 



restricts to an embedding 



Sym*{m) D' 



(m) , y ^ Yq 



(115) 
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between disks of the same dimension, which maps the geometric center —Id to 
the origin, i.e., the center of D^{m). 

Thus, by a translation and homothety inside the Euchdean space Sym{m) 
the convex subset Sym*{m) becomes, somehow, an "inward" equivariant de- 
formation of D^{m), which by (|73)) fixes the subset Sym*{m)i ~ MP"^"-*^ lying 
on the boundary sphere dD^{m). In the simplest case m = 2 we have more 
explicitly 

Sym*{2) = ^^(^l ; {u -\f + ^ (116) 

and the embedding (|115j) is actually a diffeomorphism. 

However, the embedding H115|) is not surjective for m > 2, so let us explain 
how to further deform equivariantly to make the embedding fill the whole unit 
disk. First of all, by the convexity of Sym*{m) it follows that each ray in 
D^{m) from the origin passes through a unique point of the embedded sphere 
dSym* (m). In particular, by an additional scaling we obtaing a 1-1 correspon- 
dence between the boundary spheres of the two disks, and the final composition 



: dSym*{m) ^ dD^{m) = 5"^(™+i)/2-2 ^ y ^ Yq ^ -^Yq (117) 

is certainly an 0(m)-equivariant and smooth homeomorphism. 

On the other hand, the above map may be extended to the whole disk, 
as follows. First of all, each ray from the origin intersects the embedded disk 
Sym*{m) in a segment. So, let us stretch the segment outward along the ray 
until it has unit length, but with no stretching in a neigborhood of the origin. 
Moreover, the stretching must be specified by a function on the orbit space 
since 0(m)-equivalent segments must be stretched in the same way to make V'm 
equivariant. Following this procedure we certainly obtain an 0(m)-equivariant 
and smooth homeomorphism 

: Sym*{m) L>°(m) (118) 

which extends the map in (|117() . Such an equivariant and smooth homeo- 
morhism is not unique, of course. The ambiguity lies in the group Dif fo(m) (D^i'm)) 
consisting of all equivariant diffeomorphisms of the linear model. Yet, another 
construction of equivariant homeomorphisms like V'm is decribed briefly at the 
end of this subsection. 

In summary, we arrive at the following result : 

Theorem 3.13 There is an 0{m)- equivariant and dijjerentiable homeomor- 
phism tp^ from the m-universal shape space M*{m,m) ~ Sym*{m) to its 
linear model D^{m) ~ j^m{m+i)/2~i ^ ^ _ 2 iJi^g map is a diffeomorphism. 

By referring to the following diagram 

M^[m, m) 

i TT (119) 

M < — ^ Sym*{m) — >'f'm D^{m) 
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let us explain the differential structure of Sym* (m) induced via the orbit map vr. 
Namely, the function g is said to be smooth if the composed map gojr is smooth. 
Moreover, smoothness of ip^ means the composed map o vr is smooth, and 
therefore we actually know (by our construction) that is smooth. 

On the other hand, for m > 2 there also exists a smooth map g such 
that g o is not smooth, as a function on the Euclidean disk D'^{m), and 
consequently cannot be smooth. Here we also refer to the discussion at the 
end of Section 3.4.4, together with Lemma 13.91 and Remark 13.101 Briefly, the 
differential structure on Sym* (m) induced via vr is a strict " refinement" of the 
standard structure which Sym*{m) inherits from the Euclidean disk (via ipm)- 

Remark 3.14 We have seen that the m-universal shape space Sym*(m) has 
two naturally induced smooth structures, namely induced via vr and respec- 
tively, and they are different when m > 2. In fact, Sym*{m) is not a smooth 
manifold in the first case. But in both cases vr and ip^ are smooth maps, but 
the standard structure (induced via ip^) is the only one that makes ipm 
0{m)- equivariant diffeomorphism. 

Finally, let us turn to the natural simplicial structures of the orbit spaces 
of Sym*(m) and the disk D^{m), namely the following two spherical simplices 



m—l 



Despite the above theorem, which implies that ip^ induces an 0(m)-orbit type 
strata preserving homeomorphism : A"*~^ — > A™""^, the two structures in 
H120|) are conspicuously different. The reason is that the simplicial structure 
in the first case of reflects the induced subrank stratification of Sym*{m) (cf. 
Section 3.3.3), namely the common refinement of the rank and 0(m)-orbit 
type stratification, whereas the simplicial structure in the second case is merely 
reflecting the pure 0(m)-orbital stratification of the orthogonal transformation 
group 0{m) on D^{rn). Thus, by passing from Sym*{m) to its linear model 
D^{m) the notion of rank is seemingly lost. Therefore, let us also investigate 
how the rank strata can be recognized in the linear model itself. 

Recall that M*{k, m) C Sym*[m) is the subspace lying above A'^^"'^ : r^+i = 
... = = 0, and the rank k stratum is defined by the subset Aj^~^ C A'^"^ : 
rk > 0. This stratum has all the 0(m)-orbit types 

0{ki,k2,..,kp) X 0{m- k), '^ki = k,cL^ (121) 

labelled by the various subranks k = (ki,..,kp) which record the strings of 
equalities among the numbers > 0, see (|H7j) . 

On the other hand, in the linear disk model any pattern of equality strings 
among the Xfs corresponds in the same way to a tuple k' = (ki, ..,kp,kp^i) 
of positive integers, where ki is the number of Aj's in H120() equal to Ai, k2 is 
the number of Aj's equal to A^^+i etc., and the factor 0{m — k) in (|121|) is 



39 



replaced by 0{kp^i). In any case, the last factor 0{m — k) (resp. 0{kp^i)) of 
the isotropy group H121|) is no more special than the other factors 0{ki), and 
this clearly explains why the rank is not determined by the 0(m)-orbit type. 

We claim, however, the rank k is determined in the linear model by the 
identity kpj^i = m — k provided the obvious condition ^1 ~ ^ holds. The 
latter condition merely says M*{k, m) embeds as a subset of the sphere dD^{m) 
and hence away from the interior D^{m) of the disk. The interior is, of course, 
given by ^ < 1 and here the rank \s k = m. Moreover, by removing fep+i 
from the tuple k' we are left with the correct subrank tuple k. The above claim 
about k will be settled below. 

Any 0(m)-equivariant homeomorphism 

if) : Sym*{m) D^{m) 

such as V'm; fo'^ example, induces an 0(m)-orbital strata preserving homeomor- 
phism : A™"^ A™"^. Conversely, let us see how to start from and 
construct an appropriate lifting ip as above. The idea is to construct as a 
map between fundamental domains with the 0{m)-isovariant property, that is, 
a point and its image point have the same isotropy group. Then ip will be the 
unique equivariant extension to all of Sym*{m). To choose appropriate funda- 
mental domains, first observe that A™'"^ C D^{m) is a fundamental domain 
if the numbers A, in H12U|1 are regarded as the entries of a diagonal matrix. 
Similarly, the following subset of diagonal matrices 

A™-i C Sym*{m) : Ai > A2 > ... > Ar„ > 0, ^ A^ = 1 

is a fundamental domain in Sym*{m). Then the 0(m)-orbit map projection 
A™-i A^-^ is just Ai ^/Al = ri, see JHSI, and we note that dA"^''^ 
~ ^"1-2 ^]^g subset with Am = 0. 

For example, the above equivariant homeomorphism (|118j) between the 
boundary spheres corresponds to the 0(m)-isovariant map 

cm— 2 _ _ cm— 2 

~ A'"-2 = aA™-^ (9A""-^ ~ (122) 

which sends Y = diag{Xi, .., Am-i, 0) to the matrix 

-^lo = \diag{Xi - 1/m, .., Am_i - 1/m, -1/m) (123) 
1^0 1 ^ 



where A = y Af + ... + X^-i — ^/m. 

The map ()122() is isovariant since the inequality pattern among the entries 
of the matrix Y and its image is preserved. From this it is also clear that 
Afc_|_i = if and only if the m — k last entries in H123|) are identical. This also 
settles the above claim concerning the rank k recognition in the linear model. 

3.4.4 The simplest shape spaces 

We return to the shape spaces Af*(3,m) together with their 50(3)-version for 
m = 2,3,4. Recall from Section 1 that M*(3, n — 1) = M* is the shape space 
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for the n-body problem, and for n = 3, 4, 5 the topological classification of these 
spaces are as follows : 



We shall make some further comments on the cases n = 3, 4. 
n = 3 : With the induced metric the disk is a hemisphere of the base space 

2D' = ^^^=CP' = SHl/2) (125) 

of the Hopf fibration ()112|) . which is a round sphere of radius 1/2. Its 
equator circle M*(l,2) = MP^ = 5^(1/2) represents the shapes of de- 
generate (i.e. collinear) 3-configurations. But in the study of the 3-body 
problem it is, in fact, natural to use the whole sphere as the shape space. 
The reason is that a non-degenerate 3-configuration in 3-space is geo- 
metrically a triangle which can be oriented in two different ways, by the 
ordering of the vertices. Then the two hemispheres in H125|) represent 
triangles with opposite orientation, cf. ^U], In this way a 3-body 

motion corresponds to a continuous motion of a mass triangle whose ori- 
entation changes when the motion passes through an eclipse, that is, when 
the shape curve crosses the equator circle. 

n = 4 : The boundary of the disk in (|124|) is the shape space of coplanar 
4-configurations, that is, the sphere 

M*(2,3) =CPV^2 ~ 5^ (cf. (CMl) with m = 3) (126) 

and the action of its isometry group 0(3) is equivalent (by Theorem l3.13j) 
to that of its linear model (0(3), S^, S'^p^ — !)• Namely, the linear model 
is the space of symmetric 3x3-matrices with zero trace and unit norm, 
with the natural action of 0(3) by conjugation. 

The fact that the quotient space CP^/Z2 is homeomorphic to 5"^ was already 
known to L.S. Pontryagin in the 1930's, according to ArnoldlHl, and we refer 
to [2], ^S], lU foi' different proofs of this specific result. Massey also observed 
that the induced 5'0(3)-action on CP^/Z2 has the same orbit structure as that 
of the above linear model (50(3),5^), and the existence and construction of 
an equivariant homeomorphism (such as (|117|) ) was, in fact, a problem stated 
by Massey. Moreover, Arnold[2! has constructions which are very close to our 
linear model construction. 
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Arnold also discusses the differentiable structure of 5^, as a quotient space 
CPV^2, and it is stated that CPV^2 is diffeomorphic to (cf. §1 in [H). 
However, what is shown is that the composite map 

M\2, 3) = 5^ ^ ^ CPV^2 (127) 

from to 5^ is differentiable (in the usual sense). Both CP^ and CP^/Z2 
have the induced smooth functional structure as quotient spaces of S^, but 
only CP^ becomes a differentiable manifold in this way. So, Arnold "relaxes" 
the differential structure on CP^/Z2 so that the last map in H127() becomes a 
diffeomorphism, that is, he defines the differential structure to be the standard 
structure mentioned in Remark 13.141 This makes the last map in H127|) a 
diffeomorphism while the composed map in (TTSTj) is stih a smooth 

map. The same applies, of course, to all the constructions — > S^/K — >~ S'^ 
in Section 3.4.2, where none of the "quotient spheres" S'^/K is really a smooth 
manifold when g > 4. 

Finally, let us explain why the "quotient spheres" are not smooth manifolds. 
Since all cases are analogous we consider again the simplest case 

iVr(2, 3) = S^/0{2) = CPV^2 ~ (128) 

viewed as an 0(2)-orbit space with the induced smooth structure. Then our 
claim is that the orbit space in (|128() is not a differentiable manifold in a neigh- 
borhood of the subset A = M*(l,3) ~ MP^. To see this we shall apply the 
slice theorem (cf. Section 3.1.1), according to which each 0(2)-orbit in 
belongs to a tubular neighborhood U whose image U = U/0{2) in the orbit 
space is an open set of type (|57j) . namely diffeomorphic to the orbit space of 
the slice representation. The set A represents those orbits of type 0(1), and 
the slice representation of 0(1) acts on with the eigenvalues (— 1,— 1,1,1), 
consequently any point on A has an open neighborhood diffeomorphic to 

_ 0(2) xom (M2eR2) IR2^m2 j^2 

^ ~ 0(2) - 0(1) - 0(1) ^ 



where the first factor is of type (|in2|) and is transversal to the set A and the 
second is ~ f7 H A. However, the product space in H129jl is homeomorphic to 
M^, but in view of Lemma 13.91 it is not a differentiable manifold. 



4 Geometric invariants of n-body systems 

By a geometric invariant on the centered configuration space 

M„ : ^ rmai = 

we mean a polynomial function P(ai, ...,a„) which is invariant under congru- 
ence and internal symmetries. We will describe the ring of all these invariants by 
calculating the ring of invariants for matrix spaces M{d, m) in general. Denote 
a typical matrix by 

X = (xi, ...,Xm) = [xi, ...,xj] 
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where the Xj and x* are the column and row vectors respectively. The product 
group GL{d) X GL{m) acts on M{d,m) by matrix multiplication 

X ^X^\ if) G GL{d) X GL{m) , 

and 0(d) x 0{m) is the subgroup leaving invariant the standard metric form 

m d 

I = = trace{X^X) = ^ |xip = ^ |x*|^ , (130) 

i=i j=i 

On the other hand, for fixed mass distribution /i = (mi, ..,mm), mi > 0, let us 
also consider the mass dependent metric defined by 

m 

l(i') = trace{DX^XD) = |xip , (131) 

i=l 

where D = diag{^Jmi, ^m2, .., y/rrim) and the corresponding isometry subgroup 
is 0{d) X 0*(m) C GL{d) x GL(m), where 

0*(m) = DO(m)L»"^ cGL{m) (132) 

is the subgroup leaving invariant the metric form ^ mixj on R"^ (space of row 
vectors with right side action of GL[m)). 



Lemma 4.1 The invariant ring on M = M(d,m), d <m, under the action of 
0(d) X 0*{m), is the polynomial ring with generators 

Ik= ^ mj^mj2...?njj^. |xii A Xj2 A ... A XjJ^ , I < k < d, (133) 
h<i2....<ik 

where Ii = J*^^-' and the exterior product space A^M."^ has the standard norm. 
Proof. First assume rrij = 1 for all i. It is easy to verify the identity 
Ik= Yl |xn Axi, A... Axij2 = ^ |x*^ A x*^ A ... A x*j' (134) 

ii<i2.---<ifc ii<i2....<jfc 

for each k = 1, ..,d. Then, from the column vector version it is clear that 1^ is 
invariant under the action of 0(d), and similarly the row vector expression is 
invariant under the action of 0(m). Hence, is an invariant of G. 

Let us apply the reduction principle for orthogonal transformation groups 
(cf. e.g. HHI)- Namely, we first calculate the reduced group of the action of 
G = 0(d) X 0(m) on M, which is the quotient group 

G = Ng(H)/H ~ Bd 

and is finite in our case, where 

H = AO(l)'^ X 0(m -d)cG, cf. ^ 
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is the principal isotropy group and Ng(H) is its normalizer in G. Next we 
determine the fixed point set of H 

R'^ = F{H,M) : X = diag{ri,r2,..,ra), cf. (jHU) 

which contains the fundamental domain 

M+{d) CR'^ C M :ri>r2> .... > > 

of the G-action, identified with the orbit space of G in Section 3.3.1. Moreover, 
for the G-orbit through X, rf = Aj are the eigenvalues of X^X. 

Finally, we consider the induced action of G on M*^, which is a group gen- 
erated by reflections and as such it is the Weyl group Bd, having the above 
cone M^{d) as a fundamental domain. Now, the reduction principle says that 
the restriction of polynomials induces an isomorphism between the invariant 
rings of (G, M) and (G, M'^). But the latter ring is certainly generated by the 
elementary symmetric functions 

-^1 = -^i ' -^2 = '^j'^ji -^3 = Yl ^Aj^k, .... 

i<j i<j<k 

which are also the restrictions of the functions in (|134() . and this proves the 
lemma for masses nii = 1. 

Next, let M{d,m)^'^^ be the matrix space with the mass dependent metric 
(|13H) and consider the transformation 

M{d,m)^^^ M{d,m), X XD = (^mlxi, ^/m^:>^2, V^m^m) 

where the left and right hand space have metrics (|131|) and (|13n|) respectively. 
The transformation is a G-equivariant isometry, where 0{m) acts on M(d, m)*-^^ 
by first applying the isomorphism 

0{m) DO{m)D'^ = 0*(m), cf. XU2\i 

Hence, the transformation also induces an isomorphism of the corresponding in- 
variant rings. In effect, the invariant ring of M{d, m)^^^ is obtained by replacing 
each vector Xj in (|134j) by ./rriiXi. m 

In particular, the centered configuration space M„ ~ M(3, re — 1) has the 
following basic geometric invariants, in terms of the Jacobi vectors and classical 
vector operations : 

h = J2\^if , h = J2\^i X ^jf , h = E |(xi X xj) • xfcl^ (135) 

i<j i<j<k 

for all n > 4. (For n = 3 the generators are Ii and I2 )■ The expressions for 
Ik hold for any choice of Jacobi vectors, but still it remains to express them as 
functions of ai,...,a„. To this end, consider the free n-configuration space or 
matrix space 

M„ = M(3,n)('^) = {X = (ai,a2,...,a„)} 
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with the metric form I^^'^ = |ajp as in 1)1311 By the above lemma its 

invariant ring with respect to the group 0(3) x 0*(n) has the basic generators 

Ii = J2m\^i\^ , /2 = ^ mjmj |aj X ajp , (136) 

i<j 

h= mimjmk\{a.i x a^) • a^l^ 

i<j<k 

and these are stiU complete and independent as invariants of the subgroup 
0(3) X 0*(n — 1) acting on the subspace M„ C M„, as long as 3 < n — 1. Here 
0*{n — 1) C 0*(n) is the internal symmetry group of M„ and is the subgroup 
leaving M„ invariant, or equivalently, the subgroup of 0*(n) which under the 
action on M"' (row vectors) fixes the vector (mi, m2, mn). This proves the 
following result. 

Theorem 4.2 The ring of geometric invariants of the centered n- configuration 
space Mn is the polynomial ring generated by Ii, ..,Ig in iliA'6)) . q = min{3,n — 

1}- 

Remark 4.3 Starting from the Jacobi vector expressions Ucl5\) of the invari- 
ants, a Jacobi transformation may transform them to expressions involving the 
vectors a^. In this way, however, we may arrive at a non-symmetric expression 
since ^ mjaj = 0. Below we shall give examples to illustrate the non-uniqueness 
of the symmetrization procedure and also give a geometric interpretation of In-i 
for n = 3 or 4. 

• n = 3 : Let A be the area of the triangle spanned by the vector triple 
(ai,a2,a3) with ^mjaj = 0. By simple trigonometry 

m , m fh , 

2A = — ai X a2 = — a2 x as = — as x ai 
m3 mi m2 

On the other hand, from the Jacobi vector formula 1)431) . 

xi X X2 = CiC2ai X a2 

and consequently 



/mim2m3 

xi X X2 = C1C2 ai X ^2 = 2W = A 

" m 



and 



r I |2 - fniTTlj . .2 

i2 = |xi X X2I = m |aj X aj| 

ruk 

_ mim2 I |2 , - "1-2^3 |2 , - m^rni 2 

= m |ai X a2| + m |a2 x asl + m |a3 x ai| 

3m3 3mi 3m2 

where {i,j,k} = {1,2,3} and the last sum is a symmetrization with all 
three terms equal, in fact. However, the above expression for I2 and the 
symmetric expression in (|136j) are identical as a function on M3. 
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n = 4 : Let V be the volume of the tetrahedron spanned by the vector 
quadruple (ai, a2, aa, 34) satisfying ^mjaj = 0. Then ib6V equals the 
triple product of a, — ai, z = 2, 3, 4, and 



TTi in 
6V = — |(ai X a.2) ■ asl = — |(a2 x as) • a4| etc. 

7714 mi 

Again, by expressing the Jacobi vectors Xj in terms of the a^ it follows 



xi X X2 • X3 = C1C2C3 ai X aa • ag = 3!W = V 

V m 

and consequently 
Is = (xiXX2)-X3 =m \{a.i X aj) ■ a.k\ 

m 

_ 171117121713 2 , -m27n3m4 |2 , ^ 

= m — (ai X aa) • as + m — (a2 x as) • a4 + etc. 

4m4 4m 1 

where {i,j,k,l} = {1,2,3,4} and the four terms in the last sum are 
equal. Again, Is in (|136|1 is an another symmetric expression for the same 
function on M4. 

Remark 4.4 The cases n = 3 or 4 are special since, for example, their shape 
space M* is topologically a disk or respectively. Moreover, in these cases 
M* has a unique geometric center, and this may he characterized in several 
ways. It is the unique fixed point of the symmetry group 0{n — 1), and the point 
has the same distance to each of the binary collision varieties. Furthermore, 
it is also the maximum point of I2 or respectively. That is, for a fixed size 
Ii = 1 the triangle (resp. tetrahedron) with the largest area (resp. volume) is at 
the geometric center. In terms of Jacobi vectors the condition is that Xj -xj = 
and |xj| = |xj| for i ^ j. 



5 The weighted root system of an n-body system 

A typical potential function f7(ai, a2, a„) of an n-body problem depends only 
on the pairwise distances 

rij = |aj - a^l , I < i < j < n, 

and for certain applications one would like to express them in terms of Jacobi 
vectors Xj, say 

X = (xi,..,x„_i) G Af(3,n- 1) 

is the matrix associated with a given Jacobi transformation ^ : M„ ~ M(3, n — 
1). For this purpose we shall define and investigate the weighted root system, 
depending on the mass distribution, which is primarily a book-keeping for ex- 
panding the vectors aj — a^ as a linear combination of Jacobi vectors. It is 
naturally a "weighted" version of the root system of type An-i in classical Lie 
theory, and they are identical when the mass distribution is uniform. 
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5.1 Distance functions and the ^-root system 

The mutual distances Vij have a nice geometric interpretation as functions on 
Mn which measure the distance from specific subvarieties. To explain this, 
consider the linear transformation 

M = M(3,n- 1) ^ 



defined by 



n-l 



[u,X] = ^CiXi =Xu* 



(137) 



i=l 



where u = (ci,C2, ...,c„_i) G M""^ is a given row vector of unit length and u* 
is the column vctor. Define the codimension 3 subspace B = C M to he 
the kernel 

B : [u,X] = 0. 

It is easy to verify that the distance in M from a point X to the subspace B is 
given by the function 



5b{X) = dist{X,B) = \ [u,X]\ 



*^2 
k) 



k=l 



where the are the row vectors of X and the inner product on 
standard one. 

In particular, the binary collision varieties 

Bij C Mn:ai = a^, i / j 



on— 1 



is the 



(138) 



have associated distance functions Sij of this kind, defined on M via a chosen 
Jacobi transformation ^ : M„ — > M. Consequently there is a unit vector 
(unique up to sign) 



„ C„l. „2 5 

depending on and the mass distribution, so that 

1 

3 



(139) 



*\2 
k) 



k=l 



(140) 



Let us express Sij in terms of the n-configuration X — (ai,a2, ...,a„). For 

X fixed, let Yq = (bi, b2, ....b„) be the critical point of the distance function 



5{Y) 



X| defined on B^j, where by definition bj = hj. The condition 



VS = implies 



hi = hj = 



uiia-i + mj&j f 1 _L ■ ■ 

rrii + m-i 
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and consequently in terms of reduced masses (|'29() 

S'^j{X.) = |Yo - X|^ = rrii |bj - aj|^ + mj \hj - a^l^ 
_ Tnirrij , ,2 _ 2 

mj + rrij 

that is, 

ri, = -^<5i, =-^IK,X]| (141) 
Therefore, we may define vectors Uij = —uji uniquely by the constraint 

a.i - aj = — ^ [uij,X] , i / j, (142) 

and clearly they satisfy all identities of type 

111 , , 

-Uij + ——Ujk = ——Uik . (143) 



Definition 5.1 The '^-root system with mass distribution (mi^...,m„) is the 

above collection of (2) vector pairs ib--^p=Ujj G M"^^,l < i < j < n. The 
collection {ztuj^} of unit vectors is the normalized ^-root system. 



5.2 The standard weighted root system 

It will be useful to have explicit formulas for the ^o-roots, where 

^0 : M„ ^ M(3, re - 1) 

is the standard Jacobi transformation constructed in Section 2.3. The notation 
used for various constants in Section 2.3 is also used below. We shall also use 
the notation foBy definition of ^'o (cf. (021), (jT7|) ) 

'^i ^ I \ " Sfe J- \ ^ k ■ ^ ■ 



(where the term ^Xj is undefined when j = n) and the vector Ujj = (cj^-, c^""*^) 

have its nonzero components c^- for k in the range i < k < min{j, n — 1} . In 
particular, 

U„-l,n = (0,0,...,1) (144) 

and for i < re — 1 



1 / 7re(*"i) / re^(*+i) 

:u,,,+i= (0,...,0,W j^,-J ^,0,...,0) (145) 



1 / m(* ^) I m-j+i / m„_i 

— Ui„ = (0, ..,0, 4/ rrr, J , W —7-^-7— tt) 146) 
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whereas for 2<i + l<j<n 



1 / m(*^^) I mj_|_i / mjfc 



(i) 



"^i-i / — ,0,..,0) 



For a fixed mass distribution, all ^'-root systems are, in fact, orthogonally 
equivalent to the above ^'o-root system |ib-^p=Ujj, i < j ji see the following 
subsection. 

Example 5.2 n = 3 : 



, 1712771 / mim3 

^' (m2 + m3)(mi + m2) ' V (m2 + ?n,3)(mi + m2) 



, / mim2 

ui3 = d 



(m2 + m3)(mi + m2) ' y (m2 + m3)(mi + m2) 



Example 5.3 n > 3 : 

1 



n> 3 



n > 4 



n > 5 



1 

1 



:Ui2 



m 



U23 = (0, 



U34 = (0,0, 



mi{m2 + ..+ nin) 



1712 + .. + nin 



1713 + 



rrir. 



m2[rn2 + ■■ + mn, 
m4 + .. + m„ 



■,0,...) 



m2{m3 + .. + rUn) \J 1713(1713 + .. + nin 



-,0,...) 



m3 + .. + rUn 
1713(1714 + ..+ nin) 



ms + .. + mn 
m4(m4 + .. + m^) 



,0,...) 



5.3 Weighted root systems and their metric invariants 

It follows from (|143j) - (|147|) that the (normalized) ^'o-roots Uj^ and u^; are mu- 
tually perpendicular except when {i, k} n {j, 1} 7^ 0, in which case we define 
a^''' = a'^'* to be the angle between Uj^ and u^j, namely 

cosa,'^ = Ujfc • Ufcj = -. /- ^ — — ■ (148) 

\J [mi+mk)[mk+mj) 

In particular, for equal masses the possible angles between any two non-collinear 
vectors are a = ^vr, ^vr and |7r, and the \I'o-roots constitute, indeed, a root sys- 
tem of type An-i in the usual sense. For general masses we make the following 
definition. 

Definition 5.4 A weighted root system (of type An~i) with mass distribution 
(mi, m2, mn) is a collection of nonzero vectors 

Wij = -Wji, I < i < j < n, 
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in M" ^ such that i) y/J^ij l^^ijl = 1 , ii) o-nd wji are perpendicular iff 
{i, k} n {j, 1} = 0, and Hi) the angle between Wj^ and w^j is cx^ , as defined by 

The ^'o-root system is our prototype of such a root system, and clearly the 
subset of n-1 vectors |-^p^=Uj^j_(_i| constitutes a system of simple roots for 
obvious reasons, see (|143() . (|145() . 

Theorem 5.5 For a given mass distribution, each ^-root system (see Defini- 
tion \5.1\) is a weighted root system in the sense of Definition \5.4{ and conversely, 
each weighted root system is the ^-root system for a unique Jacobi transforma- 
tion ^ : Mn M(3,n — 1). Moreover, for the given mass distribution all the 
weighted root systems are orthogonally equivalent. 

This result is a simple consequence of the fact that any ^ can be written 
uniquely as a composition 

^ = ipo^Q : Mn ^ M(3, n-1) ^ M(3, n-1) 

where if £ 0{n — 1) acts on Af (3, n — 1) by multiplication, that is, V'(^) = Xip. 
Therefore the Jacobi vector matrices X and X' of ^'o and respectively, are 
related by X' = Xip. Let {±Ujj} and |±u^^| be the normalized ^'o-roots and 
^'-roots, respectively. The defining relation (|142|) 



y^(ai - aj) = [uij,X] = [u^,X'] 

implies u^^- = Ujj^j, and consequently the \I'-root system is the orthogonally 
transformed image 

' :u'jj \ ~ \ =fc , Uij \ V (row vectors) 



I yl^i'j J I yl^i'j 

of the ^o-root system. 

On the other hand, a root system is a finite subset of M"~-^ with specified 
inner products of the vectors, and according to a result of Weyl (see Section 3.3) 
these numbers determine the subset modulo orthogonal equivalence. Moreover, 
since the weighted root system spans all M""-*^, any two of them are related by 
a unique S 0{n — 1). In particular, a weighted root system is the ^'-root 
system for a unique Jacobi transformation ^ = (p o '^q. 

5.4 On the role of the mass distribution 

The mass distribution of a n-body system manifests itself through the various 
mass dependent quantities constructed above, for example : 

• the collection {/Uj^} of reduced masses (f^ : 

• the weighted and normalized root system {itujj} ; 
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• the collection | of angles between the vectors of a weighted root 
system. 

The mass distribution (modulo scaling ) may, in fact, be reconstructed from 
any of these invariants. The reduced masses, for example, satisfy the following 
two conditions 

(0 f^i/ + f^jk > fhk ^' ^ different 

(u) Hi/ + fi^l = + fJ-Ji^ = fiii^ + f^Jk for different 

and one can define a reduced mass distribution of order n to be a collection 
of (2) positive numbers /x^^ = fij^, I < i < j < n, constrained by the above 
two conditions. Then there is a 1-1 correspondence between the usual and the 
reduced mass distributions (modulo scaling) given by the formulae for and 
their inversion 

2 

nii = ——^ — for i,j,k different 

Nj + l^ik ~ l^jk 

On the other hand, the congruence space M„ (and hence the shape space 
M*), with the kinematic Riemannian structure, is independent of the mass 
distribution. The mass distribution determines, however, the representation 
of congruence classes as points in in particular, it determines the relative 
position of the binary collision varieties. Conversely, we can reconstruct the 
masses rrii (modulo scaling) from knowledge of the position of these varieties. 

To illustrate the last statement, consider the case n = 3, where M3 is a 
round hemisphere (or sphere, see Section 3.4.4) with a distinguished equator 
circle representing the shapes of the degenerate triangles. On this circle lie 
the three collision points bi2,b23,bi3, where bij represents the shape of those 
triangles (ai,a2,a3) with = a^ 7^ and the center of mass at the origin. 
The mass distribution determines their relative position, that is, the angles (or 
distances) between the points bij, and once their position have been fixed we 
know how to determine the position of any shape (cf. jl2j). 

To describe the relative positions of the points bij, consider as the cir- 
cumscribed circle of the triangle A = A (612, 623,^31) in a Euclidean plane with 
origin at the center of S^. It turns out that A is a central triangle, in the sense 
that the center of its circumscribed circle lies in its interior. Moreover, the three 
central angles 

0</5fc<7r, ^/?fc = 2^ 

where Pj, is opposite to the vertex bij {i,j, k different), are given in terms of the 
normalized mass distribution by 



sin Pk = 27 V"^i"^2»T-3 ^ ^ 

[ruk + mi){mk + rrij) ^ 

Conversely, any central triangle can be realized in this way and the inversion 
formula is 

mk = \- ^t^^.\ (cf. CD] or [El) 
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where the positivity of m^, indeed, reflects the central property of the triangle 
A. 

The triple of angles (Pi, [^2, (3^) is in 1-1 correpondence with another triple 
(qi, a2) 03), where 

afc = 7r-/3fc/2, A; = 1,2, 3 (149) 

and these are the central angles of another central triangle, namely the tri- 
angle representing (up to congruence) the weighted root system of the mass 
distribution (cf. Section 5.3). More precisely, whereas (3)^ = P^j is the angular 
distance between bik and b^j on the circle S^, ak = o-ij is the angle between the 
normalized root vectors Uj^ and u^j since 



COSOk = Uik • Ukj = - J- ■ '—^ — ■ -, cf. (tTlHl) 

y [mi + mk)[mk + mj) 

Finally, we recall that the root system of a /c-body system is a weighted 
root system of type The way different A2 systems combine into higher 

rank systems A^, k > 3, is completely parallel to the "standard" theory of 
root systems of type Aj^., and therefore the higher rank case poses no further 
problem. 



6 Collinear central configurations revisited 

Consider the classical Newtonian potential function 



and its gradient field VU with respect to the kinematic metric © in M„ (or 
M„). An n-configuration X = (a-,^, ...,a„) is called central if VC/(X) =AX for 
some constant A, namely 

1 dU m, (ai — a,) 

Aa. = — — = ^ ' '\ i = l,...,n 151 

In fact, A = —I~'^U(yi) where / = |X|^, since by a classical result of Euler 
the homogeneity of U implies VC/(X) • X = — C/(X). The collection of central 
configurations is clearly invariant under similarity transformations. Thus we 
may fix a scaling of the vectors, say 1 = 1, and ask for solutions modulo 0(3)- 
congruence. Then the solutions are just the critical points of [/ as a function 
restricted to the shape space M*. For n = 3 the only solution which is the 
shape of a non-degenerate triangle is the equilateral triangle (by Lagrange, 
1772), whereas there are three degenerate triangle solutions (by Euler, 1767) 
and they are represented by the three Euler points on the equator circle of M| . 

Even the enumeration of all critical points in M* is, in fact, still an open 
problem for n > 4, and now the number also depends on the mass distribution. 
It is only known to be finite for n < 5. However, the number of collinear 
solutions is known to be n!/2 for all n > 3, and the first proof was presented 
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by Moulton More recently, Smale j^H] has given a topological proof using 
elementary Morse theory. In this subsection we shall give a similar and quite 
simple proof using the weighted root system of an n-body system. 

As usual, we assume that U in (|15U() is a function of n vector variables 
aj linearly related by the condition that the center of mass is at the origin. 
Actually this condition is a consequence of the identities H151|) . Anyhow, let us 
express U as a function of the n-1 Jacobi vectors 

X = (xi,X2, ...,x„_i) = [x*,X2,X3] 

defined by a fixed Jacobi transformation 

^ : M„ ^ M(3,n - 1) 

with associated root system {itwjj}. Take, for example, the standard transfor- 
mation ^'o and the expressions for the vectors Xj and Wjj in Section 2.3 and 
Section 5.2. The (row) vectors Wjj and x|. belong to the same Euclidean space 
M"^"*^ and by definition 

|aj — a^p = ^^(wjj • x^)^, Wij = Uij (cf. Section 5.1) 

k=i 

The gradient field VU is tangential to the subvariety of collinear configura- 
tions and therefore the critical points we seek are also the critical points of U re- 
stricted to the subvariety of collinear shapes, namely M*(l, n— 1) ~ RP^~^. Let 
us represent all shapes of collinear type by n-configurations X = (a^, a2, .., a„) in 
M„ with position vectors aj along the x-axis, and therefore the Jacobi matrix 
X has row vectors Xg = Xg = 0. Hence, we shall regard ?7 as a function on the 
(n -l)-space 

R"-i:x = xl = (xi,...,x„_i) 

and the condition 1 = 1 means restriction to the unit sphere S^~'^ '■ "12^1 ~ 
1. When antipodal points on this sphere are identified, we obtain the above 
projective space RP^~^, see ()154() below. 

For convenience, we shall identify the root vector Wij with the linear func- 
tional 

cJij-:M"-i^M, cJij-(x) = w.^. -x, 
and regard the family 

A = {±ujij} = A+ U A_ 
as the root system. The positive and simple roots are the subfamilies 

A+ = {uJij, i < j} D U = {uJi2,UJ23, , ,i^n--l,n} , 

respectively, with ujij = —ijJji and u}ij + ojjk = Wjfc, and following the "usual" 
procedure we define hyperplanes 

Hi^ = kerwi,- C M""^ 
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which subdivide the space into disjoint, open and connected chambers Ci whose 
union is 

Clearly, there are n! chambers, and each chamber C is distinguished by a specific 
choice of signs Sij = ±1 such that 

X gC EijUJijix.) > for all i < j 

By changing all signs one obtains the antipodal chamber. In particular, our 
fundamental chamber Cq is the one with all Sij positive, namely 

Co : iOij (x ) > for all i < j, 

or equivalently in terms of the positions aj = of the point masses (on the 
X-axis), since ujij{x ) = m — aj, 

oi > «2 > •••• > On (152) 

The above decomposition and combinatorial structure is, of course, similar to 
the Weyl chamber decomposition for the Weyl group A^-i = Sn of SU{n). 
By (|14fl)) and (|14H) . the expression (|150|) may be written 

which in the fundamental chamber reads 



i<j 



Lemma 6.1 f/(x) is a convex function on the chamber Cq, and it has a unique 
critical point xq (where U has a minimum) in the spherical chamber CqHS'""^. 

Proof. By straightforward calculations 



dU ^ wfj 

dxk f^. 'ujiji^y 

i<j 



where Wjj = {wjj,wfj,...,w^j ^) and x = (xi, j;2j j j 2;„_i). Hence, the Hessian 
of U at X S Co is the following positive definite quadratic form in the variable 

t = (^1,^2, ■■■,tn-l) 

k,i=i " ' i<j ' 

On the other hand, f/(x) > on Cq and f/(x) ^ oo as x approaches the 
walls of the chamber and is bounded away from the origin. It follows that U 
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has a unique critical point in the spherical chamber, namely a minimumspoint 
Xq. ■ 

The induced congruence group action on M"~^ is simply the group 0(1) = 
{±1} acting by inversion x — > — x, and therefore antipodal points on the sphere 
gn-2 j-gpj-gggnt collinear n-configurations with the same shape : 

In particular, a pairs of antipodal chambers is mapped to the same chamber in 
MP"~^, and consequently the latter space is divided into n\/2 chambers, each 
containing a unique critical point of C/ as a function on RP"~^. This completes 
the proof of the enumeration result originally due to Moulton. 

Remark 6.2 The above central configuration solutions (01^02, fln) can he 
distinguished by the ordering a^^ > > .... > a^^ (modulo inversion of order), 
and the minimumspoint xq in the fundamental chamber Cq is the solution with 
the ordering \15^) . However, for a given mass distribution (mi,m2, ...,mn) the 
solution satisfying Oj^ > > .... > ai„ can also be found by the same procedure, 
namely as the minimumspoint xq in the fundamental chamber Cq corresponding 
to the permuted mass distribution {rriij^^mi^, ...,mi^). In particular, for each 
string of equal masses rui = rrij = .. = the set of solutions is invariant 
under permutations of ai,aj, ..,ak. 

Let us briefly consider the explicit numerical calculation of the collinear 
central configurations by calculating the critical points of U on 5""^. Observe 
that the identity (|151j) with the additional condition 1=1 may be interpreted 
as an application of the classical Lagrange multiplier method with a constraint. 
Similarly, in the present case where we restrict to collinear configurations, we 
seek the n! solutions of the system 

V[/(x) = Ax , |x| = 1 (155) 

where A = — C/(x) (by Euler's formula) and 

The condition for a solution is the vanishing of the component of the gradient 
tangential to the sphere, namely 

VC/(x) + ;7(x)x = (156) 

As a simpleminded algorithm for finding a solution we construct a sequence of 
points Yk G S*""^ by "moving" on the sphere in the direction opposite to the 
vector (|15()|) . as follows 

_ a-U{yk))yk-^U{yk) 



\a-U{yk))yk-^U{y, 
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starting from an initial point yi in a given chamber. In general, the sequence 
will stay in the chamber and converge to the minimumspoint of U. The algorithm 
can certainly be made more effictive, but we shall not discuss these matters. 

Another approach is to introduce new variables ti linearly related to xi, a;2, ., x„_ 
via the simple roots. 

ti = uji^i+i{-x.) = Ui- Oj+i, i = 1,2, ...,n - 1 

For example, in the case n = 3 the system ()155p reads 

mi 1 — mi mi . , , 

m2 m3 

as compared to the initial system (|151|) which has three dependent equations 
similar to the above ones. In fact, the first equation of system ()15ip and (|157j) 
are identical. By elimination of A we obtain a homogeneous equation of degree 
5 which has a unique solution (modulo scaling) with both ti > 0. Thus it can be 
reduced to an inhomogeneous equation of degree 4, whereas Euler's approach 
led to the solution of a 5th order inhomogeneous equation P{uj) = (cf . pH] . 
§14). The latter equation is, in fact, obtained from the system H157() when we 
write {ti + ^2) i^t, t2 = {1 — oj)t and eliminate the resulting multiplier 

A' = \t\ 
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